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Emergence of advanced engineering materials, having remarkable technological 
characteristics in terms of high strength temperature resistance (HSTR), high hardness, 
high wear resistance, high toughness etc., plays an important role in modem manufacturing 
industries, especially, in aircraft, automobile, tool and die, medical and electronic. Shaping 
of these materials with stringent design requirements such as high precision, complex 
shapes and high surface quality is inevitable to put them in use. To achieve these 
objectives through traditional machining techniques, cutting tools with super qualities are 
needed, which are yet to be developed. Therefore, to exploit these difficult-to-machine 
advanced engineering materials with new challenges, it was perceived to develop novel 
advanced machining processes (AMPs). 

Advanced machining techniques have been classified into four basic category i.e. 
mechanical (abrasive/water jet mac hinin g (A/WJM), abrasive flow machining (AFM), 
magnetic abrasive finishing (MAP) etc.), thermal (Electrical discharge machining (EDM), 
laser beam machining (LBM), electron beam machining (EBM) etc.), chemical machining 
(CHM) and electro-chemical mac hinin g (ECM) processes. Amongst these, EDM is 
perhaps the most widely used AMP in the present day industries. EDM is a thermal 
process where material is removed by utilizing high-density thermal energy. This leaves a 
recast layer with micro-cracks on the machined surface. This continues to be the major 
concern if high surface integrity of the machined component is required. 

Technological improvement of advanced machining processes can be achieved by 
combining two or more than two different physico-chemical actions on the material being 
processed. The machining jooqqss so developed is called hybrid-machining process 
(HMP). These processes are developed to exploit the potential advantages and to restrict 
the disadvantages associated with an individual constituent process. In particular a 
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mechanical action, which is used in conventional material removal processes, can be 
combined with unconventional machining action such as mechamcal, thermal or electro- 
chemical. Electrical discharge diamond grinding (EDDG), which integrates diamond 
g rindin g and electro-discharge machining (EDM), is a new hybrid machining process for 
shaping electrically conductive very hard materials. In this process, the workpiece is 
simultaneously subjected to heating by electrical sparks and abrasion by diamond abrasive 
grains. The sparking takes place in the gap between the workpiece and metallic bonding 
material of the wheel. The improved grinding performance in EDDG is due to continuous 
in-process dressing and declogging of the wheel. In addition to this, spark discharges 
thennally soften the work material in the grinding zone and thus, dimini sh the grinding 
forces. Experimental results available in the published literature, concerning EDDG of 
electrically conducting cemented carbide, titanium alloy, and metal matrix composites, are 
technologically appealing. 

Published literature available on EDDG hybrid mac hinin g process is scarce. 
However, the available literature reveals that the feasibility of EDDG has been 
experimentally tested in cut-off and cut-off grinding configurations. The EDDG process 
can also be used in surface grinding mode with suitable design of fixture and feed control 
device. Although the investigations carried out by various researchers are valuable, but 
they seem to be rather lacking in theoretical treatment of the EDDG process. 

EDDG is a complex machining process where several disciplines of science and 
engineering are involved in its theory. It is not possible to present a simple and 
comprehensive theory, which explains the nature of the process with every detail. Analysis 
of EDDG process consists of analysis of each constituent process namely EDM and 
grinding. The work reported in this thesis makes an attempt in this direction. It comprises 
analysis of temperature distribution in the workpiece for each constituent process and also 
for combined process. Thermal residual stresses, which are induced due to steep 
temperature gradients in the worlqpiece, are the best indicators of thermal influence on the 
surfece integrity. Keeping this in view, thermal stresses durir^ cut-off grinding and EDM 
have been analyzed. The findings of this thesis will be helpful in developing a 
comprehensive surface integrity model for EDDG. 



Though simple idealized problems have exax:t solution, the real engineering 
problems like the one presented in this thesis are, in general, too complex to permit an 
exact solution- In such real problems '‘numerical approximation’ is the only recourse to 
make quantitative predictions. At present sevaal numerical mediods are available, namely, 
finite difference method (FDM), iSnite element method (FEM), finite volume method 
(FVM), spectral method, boundary element method (BEM) etc for solving different types 
of problems. Each of these techniques has its own potential and weaknesses. FEM is a 
versatile method that can solve a wide range of different engineering problems of complex 
geometry, materials, loading and boundary conditions- Ih view of these flexibilities, 
Galerkin finite element method has been employed for solving the present problem, and 
codes for temperature distribution and thermal stresses in the workpiece have been 
developed. 

In the development of thermal model for EDDG, care has been taken to make sure 
that EDM is compatible with grinding. Adjustable parameters of EDDG process are wheel 
feed rates and spark discharge parameters such as pulse on-time, duty factor, mean current 
and working voltage. Currently available thermal models of conventional grinding process 
while using conventional wheels (usually, electrically non-conducting) caimot be extended 
to grinding with diamond wheels, since behaviour of these two types of wheels is entirely 
different. With this in view, a thermal model of diamond cut-off grinding is presented to 
calculate temperature distribution in the workpiece. The model has been validated by 
theoretical /experimental data available in the literature. 

In grinding operation almost all grinding energy is converted into heat wi thin a 
small grinding zone. A large part of this heat flows into the workpiece, which results in 
elevated temperature and thermal damage to the workpiece. Keeping this in view, a 
thermal model has been developed for cut-off grinding in terms of heat conduction 
equation with specified boundary conditions. Transient temperature (primary variable) 
distribution in the workpiece is obtained by using FEM and FDM for mar ching in time. 
The depthwise variation of temperature has been compared with the results available in the 
literature for cut-off operation. Effect of feed on temperature distribution has also been 
studied. 
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Thermal analysis of EDM process has also been done by using transient heat 
conduction equation. Axisymmetric domain with relevant boundary conditions has been 
considered for analysis. In thermal modeling of EDM, spark on-/off-time, spark radius, 
heat flux profile and energy partition of discharge power are considered as the main 
parameters. The latent heat required for phase change is accounted for by modifying the 
expression for the specific heat of workpiece. Gaussian heat distribution within a spark has 
been assumed. 

Temperature distribution in the workpiece due to EDDG has been computed by 
superposition of the two temperature distributions due to grinding and EDM. In 
determining the combined total temperature, proper compatibility between EDM and 
grinding has been maintained. It has been assumed that only one spark takes place at a timft 
and next spark at the same location is repeated after sparking at all other assumed spark 
locations is over. The grinding time required to move down the grinding wheel by exactly 
one element has been completely matched with total time required for all sparks to take 
place. Results showing the effect of grindmg as well as EDM parameters on temperature 
distribution have also been reported. 

Heat generated due to grinding and EDM both will introduce thermal stresses in 
the workpiece. If these stresses are lower than the material yield stress, no thprmal residual 
stresses will be introduced. Hence, an attempt has been made to find thermal stresses due 
to cut-off grinding as well as EDM. The transient temperature distribution obtained by 
FEM is used to analyze thermal stresses in the workpiece. Plane strain condition is 
assumed for thermal stress analysis of cut-off grinding. Thermal stresses so obtained are 
below the yield strength of the workpiece material. Thermal stresses due to single spark 
have been modeled as axisymmetric problem. The zones, where thennal stresses are 
beyond yield strength are also investigated. The effects of EDM process parameters on 
thermal stresses have also been reported. 

In EDDG, it is expedient to finish-grind flat components in surface grinding mode 
with the current switched off, since the spark discharges deteriomte the surfece integrity of 
the ground surfece. In this hght, a thermal model for Electro-discharge diamond surface 
grinding (EDDSG) with zero current in which whole workpiece and a part of the grinding 
wheel are dipped in dielectric is presented by considering a moving heat source over the 
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top surface of the workpiece. The temperature field within the high speed steel woriqjiece 
has been numerically predicted for metal bonded diamond wheel by using FEM. Special 
attention has been paid to study the effects of heat flux magnitude and its distribution 
within contact length, table speed, and depth of cut on top surface temperature distribution. 

On the basis of the analysis reported in this thesis, a deep understanding about 
EDDG process can be made. The models presented are able to capture the major trends of 
surface characteristics of the machined workpiece. The theoretical approach applied is 
somewhat preliminary. As a future work, improvement can be made for the comprehensive 
surface integrity model. 

At the end it can be said that in this thesis, 3-D temperature distribution in the 
workpiece, during a new hybrid machining process called electro-discharge diamond 
grinding (EDDG), has been obtained by using 2-D temperature analysis of cut-off grinding 
and axisymmetric temperature analysis of electro-discharge machining (EDM). Thermal 
stress analysis has also been done for EDM and cut-off grinding separately. EDDSG with 
zero current has also been analyzed for temperature distribution in the workpiece. 
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Introduction and Literature Review 


1.1 Machining of Advanced Engineering Materials 

Advanced engineering materials are having greatly improved thermal, chemical and 
mechanical properties such as improved strength, heat resistance, wear resistance, and 
corrosion resistance [1]. Applications of these materials lie in those areas where high 
specific characteristics are needed. Titanium-based or nickel based superalloys, 
engineering ceramics and composites are widely used advanced engineering materials. 
Use of these materials has been escalating rapidly during the past few decades. 

Cemented carbides are enjoying an ever-increasing popularity in industry and 
especially used for metal-forming dies, cutting tools, and components requiring high 
resistance to wear. Poly cry staline diamond (PCD) is widely used for cutting tools, rock- 
drilling bits, wire-drawing dies, spray nozzle blanks, and different wear-resistant parts 
due to its unique properties of high strength and high toughness. Titanium and its alloys 
are very attractive materials for several engineering applications due to their high strength 
to weight ratio and excellent corrosion resistance [2]. The major application of these 
materials are in the aerospace industry, where it is used both in airframes and engine 
components. 

Ceramic materials are employed on account of their special thermodynamic and 
mechanical properties, their light weight characteristics, high flexural strength in 
conjunction with low density, high hardness and brittleness, accompamed by certain 
degree of chemical resistance and high compressive strength even at high temperatures 
[3,4]. Important ceramic materials are aluminium oxide (AI 2 O 3 ), zerconium oxide (ZrOi), 
titanium oxide (TiOi), Silicon nitride (Si 3 N 4 ) and silicon carbide (SiC). Silicon nitnde is 
one of the difficult-to-cut advanced materials having wide applications in gas turbines, 
rocket engines, melting crucibles, and for balls and rollers. Apart from these, other 
attractive ceramic materials are boron carbide (B4C), SiSiC and C/SiC. Particle 
reinforced metal matrix composites (PRMMCs) are a class of advanced materials with a 
wide potential for application in the automotive and aerospace industries [5]. In this 



1 


Introduction and Literature Review 


group of material, hardness, strength and resistance of the reinforcement are combined 
with the ductility and toughness of a matrix material. Aluminium is one of the most 
frequently used matrix materials due to its low density. Silicon carbide particle reinforced 
aluminium MMC is often used in the automotive industries. 

These materials are also called as advanced difficult-to~machine materials 
because of difficulty in machining due to improved material properties. Shaping and/or 
finishing of advanced difficult to machine materials with stringent design requirements 
such as high precision, complex shapes and high surface quality is inevitable to put them 
in use. To achieve these objectives through traditional machining techniques, cutting 
tools with super qualities are needed, which are yet to be developed. Machining of these 
super hard materials with complex shapes have open the way to develop novel, cost- 
effective, rapid material removal techniques that are no longer based on pure mechanical 
cutting. Therefore, to exploit these difficult-to-machine advanced engineering materials 
with new challenges, it was perceived to develop advanced machining processes (AMPs) 
[ 6 ]. 

These AMPs have been classified [6] into four basic categories. First category 
includes various mechanical processes such as abrasive/water jet machining (A/WJM), 
ultrasonic machining (USM), abrasive flow machining (AFM) and magnetic abrasive 
machining (MAM) etc. Second category includes thermal processes such as electrical 
discharge machining (EDM), laser beam machining (LBM), electron beam machining 
(EBM) etc. Third category includes bio-chemical and different chemical processes such 
as chemical milling (CHM), chemical blanking (CHB), photochemical machining (PCM) 
and bio-chemical machining (BCHM). Electro-chemical machining (ECM) is placed in 
fourth category. 

Electrical discharge machining is perhaps the most widely used AMP in the present 
day industries [7]. The scope of the EDM process ranges from the drilling of micro-holes 
that are smaller than a human hair to machining of 1,00,000-pound automotive dies. 
EDM can also be used to drill holes with high aspect ratio (upto 100:1). It is now 
believed to account for over 2% of all machining. The remarkable success of the 
application of EDM in machining can be appraised by considering the fact that many 
non-conducting materials, just for sake of machining, are rendered electrically conductive 
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to the level (>0.01 /ohm-cm) required for electro-discharge processing. It is achieved by 
doping non-conducting material with conductive elements. 

EDM is basically a thermal process [6-7] where material is removed by urilizing 
high-density thermal energy. In EDM, two electrically conductive electrodes, one being 
the tool of a predetermined shape, and other being the woriqpiece, are immersed in a 
liquid dielectric, such as paraffin or light oil (Fig. 1.1). Kerosene is a very common 
dielectric fluid used. EDM is accomplished with a system comprising two major 
components: a machine tool and a power supply. The workpiece is mounted on the table 
of the machine tool and the tool electrode is attached to the ram of the machine. The 
advancing of tool electrode into the workpiece is controlled automatically by the servo 
control system. The power supply produces a high frequency series of electrical spark 
discharges between the tool electrode and the workpiece through servo control system. 


Servo control 

Tool holder 
Tool 



Spark gap 
Workpiece 


D.C power source 


Fig. 1.1 Schematic diagram of EDM process 

During EDM, pulsed DC of 80-100 V and approximately SHz-SkHz fequency is 
passed through the electrodes, which are separated by a small gap (typically 0.01-0.5 
mm). Localized electrical breakdown of the dielectric occurs and sparks are generated 
preferentially at the closest locations (narrowest gap) between the two electrode surfaces 
where local electric field is highest. Each spark melts and vaporizes (i.e. erodes) a tiny 
amount of material from the surfaces of both the electrodes (tool and workpiece). The 
required discharge conditions are achieved by servo-controlled feed motion. Polarity of 
electrodes is selected in such a way that substantially more material is removed from the 
workpiece than finm the tool electrode. 
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There are several ways to apply spark energy for different types of operations 
depending on the type of the tool used. There are three main categories into which EDM 
based operations can be divided; sinking, cutting (including wire cutting) and grinding. 
As the name indicates electrical discharge sinking is most stiited to make irregularly 
shaped holes, slots and cavities without any tool wander. Travelling wire EDM removes 
material by electrical discharge erosion with a wire electrode traveling longitudinally 
through the workpiece. It is widely used to cut punches and dies, shaped pockets and 
many other industrial applications. Electro-discharge sawing operation is similar to 
mechanical band sawing where tool electrode is in the foim of a continuously moving 
band saw made of mild steel. 

Rotating disc EDM (RD-EDM) [8] uses a rotating electrically conductive wheel 
(similar in size to a standard abrasive grinding wheel) as the tool electrode. This process 
is used as an effective alternate for sharpening carbide and diamond tipped cutting tools, 
thus reducing the high cost of using diamond grinding wheels. Double rotating disc 
electrodes EDM (DRD-EDM) uses two rotating electrodes to erode a rotating workpiece. 
Different combinations of relative positions and angular velocities of tool electrodes and 
workpiece can produce various workpiece shapes. By keeping identical angular velocities 
of the two axes, high precision spirals can be generated with a diameter error less than 
0.004 mm, circular error less then 0.002 mm and surface roughness of Ra 63 run. 

Technological improvements of AMPs can be achieved by combining different 
physico-chemical action on the material being processed [6]. The machining process so 
developed is called hybrid-machining process (HMP). In particular a mechanical action, 
which is used in conventional material removal processes, can be combined with 
unconventional machining action such as mechanical, thennal, or electro-chemical. 

1.2 Hybrid Machining Processes 

Recently, hybrid-machining processes, which are combination of two or more 
machining processes, have attracted special interest in the field of machining advanced 
engineering materials [9]. These processes are developed to exploit the potential 
advantages and to restrict the disadvantages associated with an individual constituent 
process. Usually, the performance of hybrid process is better than the sum of their 
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performance with the same parameter settings. In some ot these processes, besides the 
performance from individual component processes, an additional contribution may also 
come from the interaction of the component processes [9]. Table l.i shows the major 
hybrid machining processes either used in industries or under development stage. 

HMPs are either electrical energy based (called as electrical hybrid machining 
processes) or mechanical energy based (called abrasive hybrid machining processes. 
Table 1 .2 shows the classification of hybrid machining processes based on their energy 
sources, material removal mechanism, tool used and medium in which process is 
performed- 

Electrical hybrid machining processes are conceived to overcome the major 
limitation of ECM and EDM in which tool and workpiece are required to be electrically 
conducting and also to increase the productivity. Electrochemical arc machining (ECAM) 
and Electrochemical spark machining (ECSM) are such hybrid machining processes in 
which the phenomenon of electro-chemical discharge is employed for material removal 
[10-12] from electrically conductive and nonconductive materials respectively In 
ECAM, discharge takes place between the tool and the workpiece due to breakdown of 
entrapped hydrogen gas bubbles causing erosion of both the electrodes. The productivity 
of ECAM is reported [9, 10] to be 5 to 50 times greater than the productivity using 
individual process of ECM and EDM. On the other hand, in ECSM the discharge takes 
place (in the form of desirable sparks) between the tool and the surrounding electrolyte in 
the vicinity of the electrically non-conducting material workpiece. 

These processes have been successfully applied for the machining of soda lime 
glass, borosilicate glass, quartz, glass fiber reinforced plastics, and ceramics [13] These 
processes can be used for hole drilling, die sinking and wire cutting for machining/cutting 
of composites and partially conducting ceramics [12]. ECAM/ECSM are specifically 
effective when dealing with materials with tensile strength higher than 1500 N/mm" and 
heat resistance alloys. Machining capacity of the order of 104 rnm^/min 

a.ccu,r2icy or 

0.04-0.02 mm and surface roughness of Ra = 1.25-2.5 mm are obtained [9] 

Ultrasonic assisted electro-chemical machining (UAECM) and ultrasonic assisted 

electro-discharge machining (UAEDM) are developed by combining ultrasonic 



Introduction and Literature Review 


machining (USM) with ECM or EDM respectively [9], UAEDM has been found [14-15] 
to increase 


Table 1.2 Gassification of Hybrid Machining Processes 


Major 
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MRR and grinding ratio while lowering grinding forces for hard materials like WC and 
TiBi. Hot machining has also been tried on advanced materials. Softening of matenals by 
localized heating with aid of a plasma jet [16] or a laser beam [17] facilitates substitution 
of grinding by other machining processes that are less expensive. 
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The abrasive hybrid machining processes are developed [9] by combining 
EDM, ECM, or USM with abrasive machining, and they are most commonly used in 
industries. These can be classified in three main subgroups: Electro-Discharge Abrasive 
Machining (EDAM), Electrochemical Abrasive Machining (ECAM) and Electro- 
Chemical Spark Abrasive Machining (ECSAM). Electro-chemical abrasive grinding 
(ECAG), Electro-discharge abrasive grinding (EDAG) and Electro-chemical abrasive 
honing (ECAH) use abrasive tool with metallic bond. Electro-chemical abrasive finishing 
(ECAF), Electro-discharge abrasive finishing (EDAF) and ultrasonic assisted electro- 
chemical machining (UAECM) use free abrasive grains. 

Electro-chemical abrasive grinding (ECAG) with metal bonded rotating abrasive 
tool is a combined process of mechanical grinding and ECM [9]. An increase in 
performance index of this process is due to improvement in surface layer properties and 
reduction in tool wear and energy consumption. ECAG is particularly effective for 
machining parts made of difficult to cut materials, such as sintered carbides, creep 
resisting alloys (e.g. Inconel and Nimonic), metallic composites (e.g. PCD-Co, AJ-SiC, 
AI-AI2CO3), etc [9]. In this process, about 90% of the material is removed by electro- 
chemical means and rest by abrasion [6], Better surface finish can also be obtained after 
machining non-homogeneous materials like cemented carbide. The quality of surface 
machined by this process is independent of mechanical properties of workpiece. Since the 
dimensional accuracy obtained in ECAG is not as good as in conventional grinding, it is 
normal practice to finish components with no-electro-chemical material removal (or zero 
current). 

Electro-discharge abrasive grinding (EDAG) is a hybrid machining process 
involving EDM and Grinding. The combination of metal bonded diamond grinding and 
EDM is the only one so far experimentally successful hence more specifically this hybrid 
process is called as Electro-discharge diamond grinding (EDDG). In this process, 
synergetic interactive effect of electro-discharge action and abrasion action are employed 
to increase machining performance. The electrical discharges of EDDG cause 
considerable decrease in grinding forces, and grinding wheel wear; and also effectively 
resharp the gnnding wheel. The abrasive action in this process helps to increase metal 
removal rate (MRR) and surface quality. The EDDG of advanced engineering materials 
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has been found to be feasible and advantageous. The details about the working of EDDG 
process are discussed in the next section. 

1.3 Electro-Discharge Diamond Grinding 

Figure 1.2 shows a schematic diagram of EDDG system, which is built on an 
EMS 5535-R 50 EZNC die sinking EDM with a provision to rotate the grinding wheel 
and thermocouple arrangement for measuring temperature in the workpiece. The rotation 
to the grinding wheel, about an axis parallel to the machine table, is given by using a 
grinding spindle assembly and a variable-speed drive mounted on the ram of the machine. 
While machining, the rotating wheel is fed under servo control. The metallic bonding 
material of the grinding wheel and workpiece surface physically separated by a gap, the 
magnitude of which depends as the local breakdown strength of dielectric for a particular 
gap voltage setting. 



Fig. h 2 Experimental set-up for EDDG L Metal bonded diamond grinding wheel 2. 
Workpiox 3. Variac It also shows the arrangement for temperature measurement 
during EDDG. 
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In EDDG process, mechanical grinding and electro-discharge erosion can be 
controlled by adjusting wheel feed rates and electro-discharge pulse parameters. The 
improved grinding performance in EDDG is due to continuous in-process dressing and 
declogging of the grinding wheel. Consequently, the grinding wheel can maintain its 
grinding ability without becoming dull. Since the contributions of mechanical abrasion 
and thermal erosion to process performances are adjustable, the EDDG can be used either 
in a grinding dominant state with a relatively less effect of electrical discharge to acquire 
a reduced heat-affected surface layer, or in an EDM dominant state with a relatively less 
effect of grinding to reduce machining force, or in a well balanced state between the 
grinding and EDM. Applications of this process include machining of components made 
of advanced engineering materials such as Al-SiC composites (duralcan), inconel, 
titanium alloy, tungsten carbide, high speed steel (HSS), cemented carbide, 
polycrystalline diamond etc. 

1.4 Modeling of Machining Processes 

The modeling of machining process can be defined from the manufacturing 
engineer’s point of view as; simplified mathematical formulation of a machining process 
which establishes a relation between input and output quantities in order to describe the 
dynamic as well as the static performance of the machining process. The purpose of 
modeling a machining operation is to develop predictive capability for machining 
performance well in advance and finally to achieve optimum productivity, quality, and 
cost [18]. 

There are many reasons why the mathematical model of a machining operation 
may be required. Some of the reasons are [19]; 

1 . Simulation of process 

2. Design and optimization of process 

3. Control of process. 

Analytical, Experimental, Mechanistic, Numerical, Artificial Intelligence (AI) and 
Stochastic methods are among the most commonly used ones for predictive modeling 
[10, 20]. In recent years, hybrid-modeling techniques, which are combination of some of 
these methods, are becoming more popular. 
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An analytical model [20] is deductively derived from basic physical principles. 
Relationships between the parameters, borrowed from physics, mechanics and material 
science, which have worked well in previous instances, are used. The main advantage of 
analytical models is that the results can easily be transferred to other machining 
conditions. Analytical models can be purely deterministic or augmented by probabilistic 
or stochastic considerations. A large number of unrealistic assumptions (to make problem 
simple and solvable) are the major limitation of analytical modeling. 

An experimental (or empirical) modeling, on the contrary, is developed by means 
ot measured values, which have been obtained in machining tests. In experimental 
modeling, the model parameters are determined on the basis of measured values and 
regression analysis methods. Every experimental model can be used for accurate 
description of the machining operation under the specified machining conditions only. 
However, compared to analytical models, the development of empirical models requires 
minimum efforts. Consequently, experimental models are used in all field of machining 
technology. 

Mechanistic methods of modeling are of more recent origin [21], This modeling is 
a combination of analytical and experimental modeling techniques. Artificial intelligence 
(Al) based modeling [21] is a new brand of modeling machining processes using artificial 
neural network (ANN), expert systems and Fuzzy Logic. Unlike analytical models which 
provide explicit models with deep physical understanding, ANN provide implicit models 
captured within the weight matrices of the net. They facilitate learning from prior 
experimental data but, not yet, from prior analytical insights. 

Another class of modeling methods, which is the only recourse to make 
quantitative predictions of complex real engineering problems, is numerical method 
[18,21], At present several numerical methods are available, namely, finite difference 
method (FDM), finite element method (FEM), finite volume method (FVM), spectral 
method, boundary element method (BEM) etc for solving different types of engineering 
problems. Each of these techniques has its own potential advantages and weaknesses. 
FEM is a versatile method that can solve a wide range of different engineering problems 
of complex geometry, heterogeneous materials, varying loading and different boundary 
conditions. In view of these flexibilities, FEM has been employed for modeling of 
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machining processes in recent years. The major limitations of applying FEM for 
modeling of machining processes are large computational burden and incorporation of 
dynamic aspects. 

The foregoing discussion clearly indicates that modeling of EDDG process for 
various aspects are needed to obtain deep understanding about the process behaviors. The 
following section provides a review of the relevant literature available. 

1.5 Literature Survey 

The Electro-Discharge Diamond Grinding (EDDG) has been developed in USSR 
during late seventies. But in late 80s its applications gained momentum in various 
industries where there is a need to grind parts made of difficult-to-machine materials. 
EDDG process is in its early stage of development and publications available are scarce. 
Most of the researchers have conducted experiments to investigate the effects of process 
parameters on its performance. Theoretical treatment of EDDG process is almost nil. An 
extensive literature survey yield few references related to experimental studies on EDDG 
and theoretical treatment of its constituent processes, e.g. EDM and Grinding. A brief 
review as given below, is presented in three parts, viz. experimental investigation of 
EDDG, thermal modeling of grinding and thermal modeling of EDM. 

1.5.1 Experimental Investigations of EDDG 

Grodzinskii and other researchers [23-26] have done extensive experimental work 
on EDDG But they all are more of an exploratory in nature. The role of electrical 
discharges on grinding forces, grinding wheel wear, and geometrical accuracy while 
grooving and cutting-off of cemented carbide and few other advanced ceramics have 
been experimentally studied by Aoyama and Inasaki [27]. They used electrically 
conducting diamond grinding wheel which is made by making grooves around the 
periphery of the wheel to insert resinoid bonded diamond grains of 100 jam grit size. 
They found that the normal and tangential grinding forces decrease with an increase in 
the applied voltage at the expense of an increase in wheel wear. The reasons for decrease 
in the grinding forces have been attributed to two factors (i) some portion of the uncut 
material is removed by sparks erosion, and (ii) self-sharpening of the grinding wheel due 
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to sparks. The primary limitation of the process was found to be high wear rate of the 
wheel. 

Rajurkar et al. [28] have reported the characteristics of EDM-grinding hybrid 
process, which they call Abrasive electro-discharge grinding (AEDG). Experiments have 
been conducted to machine AJ-SiC composite, and titanium alloy with copper bonded 
diamond wheel. The electrical parameters, such as peak current and pulse on-time, and 
wheel speed as mechanical parameter have been identified to be the main parameters 
influencing the process responses for a particular wheel and servo reference voltage. 

Introducing abrasion effects into the EDDG process has been observed [28] as 
responsible to increase in material removal rate in comparison to conventional EDM and 
EDG [Fig. 1.3a]. Material removal rate of the EDDG process at pulse on-time of 10 ps 
and wheel speed of 500 rpm is found to be about five times greater than that of EDM and 
about twice of EDG in case of Al-SiC composite machining. Effect of pulse on-time on 
material removal rate is not significant [Fig. 1.3b] during machining of Al-SiC 
composite. But in case of EDG and AEDG, under the same conditions, a material 
removal rate is markedly increased (Fig. 1.3b). It may be an evidence of better utilization 
of electrical discharge energy. Similarly, higher peak currant yields higher removal rate 
due to the intensified electro-erosion (Fig. 1.3c). 

Koshy et al. [29] have found a way for removing problem of rapid wear of 
diamond wheel due to graphitization during grinding of ferrous materials (HSS) by using 
electrical discharge diamond grinding (EDDG). The role of current and wheel speed on 
the material removal rate, the grinding forces and power consumption have been studied 
to elucidate the mechanism of material removal. The trend exhibited by the MRR-wheel 
speed characteristics has been found to be different in different current ranges [Fig. 1.4a]. 
They have found that MRR initially increases with wheel speed but tends to decrease at 
higher end of the wheel speed at a current of zero and lA value due to blunting of the 
abrasives. But the MRR increases monotonically with speed at the current of 2.5 and 5 A. 
They have found that specific grinding energy in EDDG decreases with increase in 
current for any wheel speed [Fig 1 .4b]. The effect of current during EDDG of cemented 
carbide has also been studied [30] [Fig. 1. 5a]. Such studies are advantageous for exploring 
the possibility of use of less expensive abrasives in place of costly diamond wheel. 
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Choudhury et al. [31] have studied the effect of current on MRR and grinding 
forces for different voltage, pulse on-time and duty factor during EDDG of HSS. It has 
been observed that tangential grinding force decreases with increase in voltage and duty 
factor for a particular value of current [Fig. 1.5b]. They have also reported the effect of 
process parameters on the MRR. 



(a) 



(b) 


(c) 


Fi^re 1.3 (a) Effect of applied voltage on specific grinding forces (b) Effect of pulse 
on time on MRR for different wheel speed and (c) Effect of pulse on-time on MRR for 
different current, in EDDG I28J 
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Figure 1.4 (a) Effect of wheel speed on MRR for different current (b) Effect of 
current on the specific grinding energy for different wheel speed during EDDG [29] 



(a) (b) 

Figure 1.5 (a) Effect of current on normal force for different pulse on-time 
(cemented carbide workpiece) (30] (b) Effect of current on tangential force for 
different duty cycle (HSS workpiece) (31J during EDDG 
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1.5,2 Thermal modeling of Grinding 

The grinding process is an important precision machining process, which requires 
high specific energy input in comparison to other machining processes. In any grinding 
process, virtually whole of the grinding energy is converted into heat in the grinding zone 
(wheel-workpiece contact) due to friction and plastic deformation. This leads to high 
temperatures in the workpiece which can cause undesirable effects at the ground surface 
such as thermal damage (bum marks and cracks), variation in hardness, and generation of 
residual stresses. Elevated temperatures generated in the grinding zone also have essential 
influence on the performance of the grinding wheel. In an attempt to avoid these effects, 
it is important to be able to predict workpiece temperature distributions in terms of 
grinding conditions, so that the process parameters can be adjusted to yield acceptable 
workpiece temperatures. Many researchers have intensively studied thermal modeling of 
grinding and many papers have been published. 

Most of the thermal models developed for grinding are mainly focused on heat 
transfer to the workpiece with specified heat flux and with or without convection at the 
workpiece surface. 1 he earlier works related to thermal model of grinding have been 
compiled in detailed reviews by Snoeys et al. [32], Malkin [33] and Zhang et al. [34]. 
Some findings related to the development of thermal modeling of grinding are listed in 
Table 1.3. The pioneering and widely accepted work for determining temperature field in 
shallow cut grinding was carried out by Jaeger [35], who presented a two-dimensional 
model with a band heat source of constant strength moving with a constant velocity 
across adiabatic semi-infinite stationary body. Most researchers up to the present day use 
this model with little modifications to make it more realistic. 

Cutwater and Shaw [36] assumed that the grinding energy is generated at the 
shear plane only. No consideration was given to the heat generated at the wear flat areas. 
Also, with cutting edges of large negative rake angles, this model under-estimate the 
proportion of the heat energy entering into the workpiece. In contrast, Hahn [37] pointed 
out that the friction forces at the wear flat areas (grain-workpiece rubbing surface) is 
extremely important and it is the principal heat generation source in grinding. The heat 
energy generated at shear plane can be neglected. 
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Table IJ Key developments in thermal modeling of Grinding 


Year 

Author(s) 

Ref. No. 

Remarks 

1942 

Jaeger 

35 

Sliding uniform rectangular heat source model. No 
cooling. 

1952 

Cutwater and Sliaw 

36 

Sliding heat source at the shear plane. No cooling 

1956 

Halm 

37 

Sliding heat source at the grain-workpiece rubbing 
surface (wear flat). No cooling 

1970 

Oes Ruisseaux and 
Zerkle 

38 

Temperature at workpiece surface which remains after 
grinding due to shear heat source. Jaeger’s nKxlei is 
applied. No cooling 

1970 

Dcs Ruisseaux and 
Zerkle 

39 

Sliding uniform rectangular heat source model with 
uniform convective cooling 

1974 

Malkin 

40 

Temperature rise due to the frictional energy at the wear 
flats. 

1978 

Snoeys et aJ. 

32 

Vanous heat flux model-triangular, rectangular and 
parabolic 

1983 

Peters and Vansevenant 

41 

Jaeger’s model supplemented with the consideration of 
heat loss to the chip, wliich is valid for creep feed and 
pendulum grinding. 

1984 

Kopaiinsky 

42 

Temperature distribution by super imposing the two 
distributions (1) produced by grits which have left the 
contact area; (2) individual gram cutting edge with wear 
flats ourently engaged with workpiece. 

1988 

Lavine 

43 

Convective heat transfer between tlie workpiece, wheel 
and grinding fluid, wheel and grinding fluid together is 
called composite solid. 

1991 

Lavine 

i 44 

Major part of the grinding energy is generated at wear 
flat area and some part is generated at the shear planes. 

1994 

Guo and M^ilkin 

45 

Tlierraal model to predict the burnout heat flux liimt at a 
critical temp, for film boiling. 

1994 

Zheng and Gao 

46 

Thermal model for slotted / segmented wheel. Heat flux 
shape- triangular, rectangular and trapezoidal 

1995 

Jen and Lavine 

47 

Grinding zone temperature and heat flux distnbutions, 
DuhmaJ’s theorem is used. No need to specify a priori 
the value of the fraction of grinding energy entoing the 
workpiece or the convective heat transfer coefficient of 
the grinding fluid 

1995 

1995 

Zhang and Mahdi 

Guo and Malkin 

48 

49 

FEM based steady state thermal model with moving 
hftat flux of triangular profile. Convective cooling. 

FDM based thermal model to find transient temperature 
distribution during cut-in. steady state and cutout 

1996 

Zhang and Faghii 

50 

Thermal model for grinding zone heat flux and 
temperature distributions with or without film boiling. 
Integral approximation method is used. 

1997 

Biermaim and Schneider 

51 

52 

Workpiece unsteady temperature dismbutioit 
rRotangiilar heat source and convective cooling 

2000 

Demetriou and Lavine 

Coupled thermal model of the wheel grains, fluid, chips, 
and woriqjiece to predict heat flux into the woikpiece 
and resulting temperature distribution. Dnhmal’s 
theorem is used. 







18 


Introduction and Literature Review 


The models [35-38] discussed above are used for predicting workpiece surface 
temperatures without considering the effect of surface cooling (dry grinding). The effect 
of surface cooling by the grinding fluid was analyzed by Des Ruisseaux and Zerkle [39], 
They have considered constant convection coefficient over the entire top surface of the 
semi-infinite body. It was found that, convective cooling outside the contact region is 
more important than within the contact region for bulk cooling of the workpiece. 

Malkin [40] assumed uniformly distributed heat flux due to frictional energy at 
the wear flats to find temperature rise over the workpiece surface in dry grinding. It 
implicitly indicated that frictional heat generated at the grain-chip interface was relatively 
small in comparison to wear flat area. Snoeys et al. [32] have argued that transient effects 
should be negligible if the workpiece speed is higher than about 50 mm/s for grinding of 
steel. This fact is generally true for regular grinding but not necessary for creep-feed 
grinding. Kopalinsky [42] has concentrated on predicting workpiece temperature 
distribution by superimposing the two temperature distributions, one produced by grits 
which have left the contact area and other produced by wear flats currently engaged with 
workpiece. 

All the papers discussed above have assumed some value for the fraction of 
energy entering the workpiece. Lavine [43] eliminated the need to specify a priori the 
fraction of the energy entering the workpiece, by developing a model of convective heat 
transfer between the workpiece, wheel and grinding fluid. She has considered wheel and 
grinding fluid together as composite solid. However, it was assumed that all the grinding 
energy was generated at the wear flats rather than at the shear planes. But, she [44] 
modified her heat transfer model to account for the grinding energy generated at the shear 
plane. With this consideration, the results of conventional grinding with CBN wheels 
were found significant but not for AI2O3 wheel. 

Guo and Malkin [45] found the bum-out heat flux limit at a critical temperature 
for film boiling during creep-feed grinding. Beyond this limit cooling becomes 
ineffective and temperature rises catastrophically. Zheng and Gao [46] also suggested use 
of slotted or segmented wheel as one of the effective ways to reduce grinding thermal 
damage. They developed a more general thermal model for prediction of grinding 
temperature for segmented or slotted wheels. 
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Jen and Lavine [47] developed a thermal model in down grinding using 
Duhamel’s theorem to solve the conjugate problem of heat transfer in the workpiece, 
wheel grains and fluid, in order to predict the workpiece surface temperature in the 
grinding zone. Zhang and Mahdi [48] used finite element method to simulate surface 
grinding process. Temperature distribution was obtained considering triangular heat flux 
profile moving with uniform velocity over the workpiece top surface. Convective cooling 
was also taken into consideration, Guo and Malkin [49] used finite difference method to 
find transient temperature distribution in the workpiece during cut-in and cut out period. 
They have also found that the workpiece actually follows a transient behavior during the 
entire grinding pass for thermally short workpieces. 

Zhang and Faghri [50] developed a heat transfer model to simulate the heat 
transfer to the workpiece, wheel, and fluid during surface grinding using integral 
approximation method. The grain is approximated as semi-infinite fhistum of cone for 
finding heat transfer to the abrasive grain. This model was also used to simulate the 
grinding process when the boiling occurs in the grinding zone, in which the heat flux has 
discontinuous changes. 

Biermann and Schneider [51] used Jaeger’s model and finite element method to 
simulate the surface grinding of cemented carbides with resin bonded diamond grinding 
wheel. Also, the numerical results were supplemented by measurements of the grinding 
temperature. Recently, Demetriou and Lavin [52] developed a completely coupled model 
of the wheel grains, fluid chips, and workpiece with which it could be possible to predict 
heat flux into the workpiece, and the resulting temperature distribution. 

Most of the thermal models developed as discussed above are applicable for 
surface grinding. References available on thermal modeling of cut-off grinding are not 
many. 

Eshghy [53] was theoretically analyzed the thermal aspect of abrasive cut-off 
operation to find the temperature of cutting zone. The mathematical model developed by 
him consists of two interacting one-dimensional semi-infinite domains m transient 
conditions. The solution was obtained by Laplace transformation to find cutting zone 
temperature. Korolev [54] studied the temperature distribution in abrasive cut-off 
operation with rubber bond wheels of 2-4 mm thickness cutting 10-mm long rods of steel 
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at peripheral speed of 50 na/s. The mean surface temperature at the wheel-work contact 
layer of metal, termed as contact temperature, was measured by using thermocouples. He 
found that the temperature rises rapidly as the wheel penetrate the first few mm into the 
steel workpiece, and drop steadily as the wheel heats up causing bond softening. Korolev 

[54] concluded that various cutting conditions that promote self-sharpening of the wheel 
lowered the contact temperature. 

1.5.3 Thermal Modeling of EDM 

Electric discharge is a complex phenomenon of EDM process where several 
disciplines of science and engineering are involved in its theory. Researchers have 
proposed several mathematical models based on heat transfer techniques. Some key 
findings in the development of thermo-mathematical modeling of EDM are listed in 
Table 1.4 and also discussed in the following section. Melting and vaporization of 
material are considered to be important phenomenon in predicting MRR in EDM process 

[55] . 

Van Dijck and Dutre [56] developed a thermal model based on two-dimensional 
transient heat conduction in a semi-infinite body subjected to a time-dependent uniform 
heat flux over a small circular part of its upper surface. The computational procedure for 
numerical calculation of crater boundary was based on the exact solution methodology. 
Toren et al. [57] modeled electrical erosion based on the propagation of the melting 
craters with evaporation. The discharge channel was considered as constant point heat 
source. Numerical quadrature rule was used for finding melting radius and evaporation 
radius. Marty [58] found that a thermal model based on a disc heat source on the surface 
of a semi-infinite solid with a constant heat flux is a good approximation for the spark 
erosion process. The crater depth was considered as a function of both the pulse duration 
and radius of heat source at work surface. 

Erden and Kaftanoglu [59] presented a thermal model based on the uniformly- 
distribution-point-source assumption. In this model, the main assumption is that the 
electrode surface intersecting the spark channel is covered by point sources. An increase 
in discharge radius means that there are more point sources over the heat source area. 
Total energy released by these sources is equal to total spark energy. It was also reported 
[59] that only 10-20% of the total molten mateiial is removed at the end of a pulse. 
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plasma model was also developed [67] for spark created by EDM in a liquid media. 
Madhu et al [68] developed a thermal model using finite element method (FEM) for 
determining material removal and depth of damaged layer during EDM. They assumed 
Gaussian heat distribution of power within a spark and updated the machined surface 
periodically. 


Table 1.4 Key developments in thermal modeling of EDM 


Year 

Author(s) 

Ref. 

Remarks 

1971 

Snoeys and Van 

Dvck 

55 

Semi-infinite cylinderical workpiece. Disc heat source 

1974 

Van Dizck and Dutre 

56 

Finite or semi-infinite body, Time dependent uniform heat flux 
over a small circular part. 

1975 

Toren et al. 

57 

Heat conduction model with constant point heat source. Melting 
radius and evaporation radius is calculated. 

1977 

Marty 

58 

Constant disc heat source over a semi-infinite solid. Crater depth 
and radii of heat source are correlated by an universal function 
for given heat input and maximum surface temperature j 

1980 

Erden and 

Kaflanogiu 

59 

Two themial models are developed based on (1) uniformly 
distributed point-heat-source (2) Variations of heat source 
densih' and its radius. Agreement with experimental results arc 
closer with second model but computational time is more. 

1980 

Erden and 

Kaflanogiu 

60 

Optuniration of discharge energy pulse form 

1980 

Dharmadiiikan and 
Sharma 

61 

Multi heat source model. 

1982 

Jilaniand Pandey 

62 

2-D heat conduction model, Seim-infinite body. Disc heat source 
with constant radius 

1986 

Pandey and Jilani 

63 

Semi-infinite body, disc heat source with expanding plasma 
channel. 

1983 

Pandit and Rajurkar 

64 

Tliermal model based on data dependent systems (stochastic 
empirical). Predictions are much closer to the experimental 
results 

1989- 

93 

DiBitonto and Co- 
workers 

65-67 

Point heat source model for (PHSM) cathode erosion. Expanding 
circular heat source model (ECHSM) for anode erosion. Also, a 
comprehensive spark model for the EDM plasma based on 
equations from fluid mechanics, thermodynamics and radiation. 
Superheating mechanism of EDM is explained. 

1991 

Madhu et al. 

68 

REM based heat conduction model with Gaussian heat 
distribution of power within a spark and periodic updating of 
machined surface. Heat absorption owing to melting is also 
accounted 

1996 

Bhattachaxya et al. 

69 

Model used m ref. [68] is updated by including more realistic 
effects such as on/-off-time and latent heat of melting and 
vaporization. 

1997 

Shankar et al. 

70 

An integrated approach where discharge phenomenon and 
temperature distribution in the electrodes have been treated 
together. Concave spark profile in the centre is obtained 
theoretically. 
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Further, Bhattacharya et al. [69] improved this model by accounting more realistic cases, 
such as on/-ofF-time and latent heat ot melting and vaporization. An integrated approach 
has also been applied [70] where the discharge phenomenon and temperature 
distributions in the electrodes have been treated together. Spark profile across the gap in 
terms of the spark radius at different cross-sections has been found to be different 
(concave at the centre). 

Survey described above show that EDDG is a complex hybrid machining process 
where several disciplines of science and engineering are involved in its theory. It is 
difficult to formulate a mathematical model, which explains the nature of the process 
with every detail. This thesis makes an attempt to develop a theoretical thermal analysis 
model of EDDG process. The scope and objectives of the present work, and the 
organization of the thesis are described in the following section. 

1.6 Scope and objectives of the present work 

Published literature available on EDDG hybrid machining process is scarce. 
However, the available literature reveals that the feasibility of EDDG has been 
experimentally tested in cut-off grinding and cut-off grinding configurations. Most 
previous studies of EDDG are experimental. Demands for theoretical studies have arisen 
because purely experimental investigation is too expensive and some times very difficult 
to perform. A theoretical method can simultaneously account for many principal factors, 
like those from mechanical, thermal and chemical effects, so that it is captured in terms 
of physics and mathematics. 

Thermal modeling of EDDG is of major importance in order to avoid thermal 
residual stresses and thermal damage to the workpiece. Hence, thermal model can be 
developed to determine temperature distribution within the workpiece. The residual 
stresses in the workpiece are the best indicators of thermal and mechanical influences on 
the surface integrity. The residual stresses are generated in the workpiece due to plastic 
deformation, thermal disturbances, and phase transformation during EDDG. Hence, 
thermal stress model is needed to determine stress distribution within the workpiece. 

The work reported in the thesis comprises analysis of temperature distribution in 
the workpiece for each constituent process of EDDG and also for combined process. 
Thermal stresses, which are induced due to each of the constituent process of EDDG, 
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have also been analyzed. Also, temperature distribution behaviour in the workpiece has 
been analyzed during EDDG when it is applied in surface grinding mode without spark to 
understand the nature of variation of temperature in diamond surface grinding. The 
findings of the thesis will be helpful in developing a comprehensive surface integrity 
model for EDDG. 

The following are the objectives of the present work: 

■ To analyze thermal aspects of cut-off grinding and EDM (both are the constituent 
process of EDDG) using finite element method. Rectangular domain for diamond 
cut-off grinding and axisyemmtric domain for EDM are considered for 
temperature analysis. 

• To develop and analyze thermal model for EDDG using superposition theorem. 

■ To develop mathematical models for determination of components of thermal 
stresses in cut-off grinding and EDM both. 

■ To analyze the thermal aspects of Electro-discharge diamond surface grinding 
(EDDSG) with zero current. 

L7 Organization of the thesis 

First chapter of the thesis comprises introduction of hybrid machining processes 
with their potential and limitations, working principle of electro-discharge diamond 
grinding, need of modeling of machining processes, potential of FEM, review of 
literature, and scope and objectives of the present work. Thermal analysis of cut-off 
grinding and EDM with associated results of temperature distribution only in the 
workpiece constitutes chapter 2 of the thesis. The results of superposition of the two 
temperature distributions to get temperature distribution in the combined process 
(EDDG) are also reported in the same chapter. Chapter 3 presents two thermal stresses 
models in which the first one is used for cut-off grinding and the second one is used for 
EDM. Thermal analysis of electro-discharge diamond surface grinding (EDDSG) with 
zero current is presented in chapter 4. Chapter 5 presents the main conclusions of this 
work and guidelines for future work. 



Chapter 2 

Thermal Analysis of EDDG 
2.1 Introduction 

Electro-discharge Diamond grinding (EDDG) integrates diamond grinding and 
EDM. Heat energy generated during EDDG (due to grinding and EDM both) can cause 
undesirable effects such as micro-cracks, generation of critical residual stresses, and 
metallurgical damage that may affect the surface integrity of finished components. An 
accurate prediction of thermal history (or transient temperature distribution) in any 
machined component is a prerequisite for a reliable prediction of surface integrity. But 
hardly any analytical/empirical/numerical study has been reported in the literature for 
predicting transient temperature distribution in the workpiece during EDDG process. 

Modeling of such a complex process like EDDG is not possible by analytical 
method because of complex boundary conditions and coupling of EDM and grinding 
processes. EDDG process is too complicated to warrant appropriate analytical models. 
FBM is a versatile method that can solve a wide range of different engineering problems 
of complex geometry, heterogeneous materials, varying loading and different boundary 
conditions. In view of these flexibilities, FEM is employed for modeling of EDDG 
process. 

Configuration of a typical EDDG process in cut-off grinding mode is shown in 
Fig. 2. 1 . Here, L is workpiece length, Vf is feed velocity, Vs is wheel speed, D is diameter 
of metal bonded diamond grinding wheel, B is width of workpiece and H is height of the 
workpiece. Two electrodes of the EDDG namely workpiece and metal bonded diamond 
grinding wheel are dipped in liquid dielectric (usually kerosene) as shown in Fig. 2.1. 
Metallic wheel bond and work surface are physically separated by a gap, magnitude of 
which depends on protrusion height of grains and uncut thickness as shown in Fig. 2.2. 

In the present chapter an attempt has been made to determine temperature 
distribution in the workpiece during EDDG for metal bonded diamond wheel in cut-off 
grinding mode (Fig. 2.1). Finite element based codes have been developed for the 
simulation of temperature distribution. Separate codes have been developed for grinding 



Fig. 2.1 Schematic of EDDG process in cut-off grinding mode 



» yrnr 




Thermal Analysis of EDDG 


27 


and EDM due to the difference in the nature of the temperature distribution, i.e being 2-D 
for grinding and axisymmetric for EDM. The temperature distribution in whole 
workpiece domain due to EDDG is obtained by superposing the temperature distributions 
obtained from these two codes. 

The workpiece is simultaneously subjected to heating due to electrical sparks 
generated between the bond and workpiece, and heat generated due to abrasion by 
abrasive grains having protrusion height more than the inter-electrode gap (Fig. 2.2). 
Both these phenomena (sparking and location of a grain and thereby abrasion), under 
actual operation, occur randomly over the top surface of the workpiece. The transient 
temperature distribution due to EDDG process is considered as superposition of transient 
temperature due to both-grinding heat source as well as EDM heat source neglecting the 
interaction effect between grinding and EDM. The interaction effect may be in terms of 
reduction in grinding forces because of temperature softening effect due to EDM and 
increase in material removal rate by EDM due to flushing action of grinding. The stresses 
developed during machining will also change if interaction effect is taken into 
consideration. Here, temperature distribution due to grinding is determined by 
considering whole domain as 2-D boundary-initial value problem whereas, temperature 
distribution due to sparking action is treated as axisymmetric boundary-initial value 
problem. 

Assumptions 

Due to random and complex nature of EDDG process, following assumptions are made to 

make the problem mathematically tractable. 

1. The workpiece material is homogeneous and isotropic. Although the materials 
used in real practice like HSS are not homogeneous in structure however to 
simplify the problem it is considered to be homogeneous throughout Hence, the 
average values of thermal properties are used. 

2. The thermal properties, viz. specific heat and thermal conductivity, of the 
workpiece material are treated as temperature independent. 

3. Only a fraction of grinding as well as spark energy is dissipated as heat into the 
workpiece. Rest of the heat is assumed to be distributed between the grinding 
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wheel, chip, debris, and dielectric. However, the present study is restricted only to 
the analysis of temperature distribution within the workpiece. 

4. The protrusion height of all the grains is equal and remains constant throughout 
the operation. Since the time of computational machining taken in this model is 
short (=0.05 s) hence grain wear, fracture and change in grain spacing for this 
small period of machining time can be considered as negligible. While 
considering machining for a longer duration, grain wear and fracture are required 
to be considered with probabilistic approach which can account the randomness 
of grains as well as their behavior. 

5. Length scales of grinding and EDM domains are different. Grinding 2-D domain 
considered in present study is 6 mm X 20 mm whereas EDM axisymmetric 
domain is considered within 2 mm diameter and 1 mm depth. 

The above stated assumptions 1-3 are valid for both the cut-off grinding as well as EDM 
where as assumption 4 is valid for cut-off grinding only. Assumption 5 is used for 
superposition. These five assumptions are not valid in practical cases. Since it is the first 
theoretical approach to model a complex hybrid process like EDDG, some assumptions 

are made to simplify the problem. Otherwise the problem will become complex which 
will be difficult to solve. 

2.2 Thermal Analysis of Cut-off Grinding 

In cut-off grinding operation, which is a constituent process of EDDG, it is 
assumed that total grinding energy is converted into heat within a small grinding zone, 
which leads to high temperature in the grinding zone. The heat is taken away by the 
workpiece, grinding wheel, dielectric, and chips. To avoid thermal damage to the 
workpiece, it is important to understand heat transfer in the workpiece. Here, only heat 
transfer to the workpiece is analyzed. Some researchers have intensively studied heat 
transfer in grinding experimentally as well as analytically (see section 1.5.2). Heat 
transfer study has been done only analytically for abrasive cut-off operation [53-54,71] to 
find the cntical down feed (easiest one to vary) for maximum temperature for different 
wheel work combinations. Shaw et al. [71] found that specific energy (U^) decreases with 
down feed in cut-off grinding operation. 
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2.2. 1 Mathematical Modeling 

In this section, a thermal model for the cut-off grinding is presented. Knowledge 
of operating temperature within the grinding wheel and the workpiece is essential to 
understand the cutting phenomenon involved. To find temperature distribution in the 
workpiece, proper modeling of heat flux profile and energy partition for workpiece- 
abrasive wheel combination are needed. 

The heat flux generated during cut-off grinding is considered as a rectangular 
heat flux over the workpiece top surface (Fig. 2.3). Here, L and H are the length and 
height of the workpiece enclosed by the boundary surfaces SI, S2, S3, and S4. Its value is 
calculated by using following equation [72]. 

a = R U Vr (2-^) 

where, Rwg is energy partition due to grinding, Vf is feed given to the workpiece, and Uo 
is specific grinding energy. 



Heat Source 


H 


X 


Fig 2 J Thermal model of cut-off grinding 
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The effect of cutting speed (V,) on qwg is taken care of by specific grinding 
energy (Uo) term because they are related with each other by following equation [72]. 

FV FV 

U =-L£-i- = -S-L (2.1a) 

" MRIi V^A 

Here, Fc is tangential grinding force, Vj is grinding speed related to RPM of grinding 
wheel, MRRv is volumetric material removal rate and A is the area of cross-section of 
workpiece top surface. 

Selection of energy partition (Rwg) value, which is defined as a ratio of the 
energy going to the workpiece and the total grinding energy, is of special importance in 
determining thermal load on the workpiece. Shaw [71] has reported that the value of Rwg 
decreases from 0.31 to 0.22 by increasing feed rate from 4.23 to 12.7 mm/s during dry 
cut-off grinding of AISI 1020 steel workpiece with alumina grinding wheel. It further 
reduces (upto 0.10) with feed if special graded alumina wheel is employed. But, 
information about Rwg for metal bonded diamond wheel in cut-off grinding of steel is 
nonexistent in the literature. Here, the value of Rwg is chosen as 0.08 for cut-off grinding 
of high speed steel (HSS) workpiece with metal bonded diamond wheels, which is closer 
to the value of Rwg for special graded alumina grinding wheel and AISI 1020 workpiece 
[71]. Further, high thermal conductivity of diamond grains and metallic bond of a 
diamond grinding wheel leads to slightly less energy partition (< 0.10) than used for 
special graded resin bonded alumina grinding wheel. 

2.2.2 Finite Element Formulation 

The physical phenomenon of cut-off grinding can be described by boundary- 
initial value problem i.e. by energy equation (which is a differential equation) and 
appropnate boundary and initial conditions. Solving this problem using available 
classical methods would be very difficult and tedious. Finite element method can be used 
for such complicated problem to obtain approximate solution by converting its governing 
equation (2.2) into set of algebraic equations. 
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Governing equation 


Temperature field due to cut-off grinding is governed by thermal diffusion 
equation; 


dT d dT ^ d dT ^ 

dt ox ox Oz oz 


( 2 . 2 ) 


where, T is temperature, t is time, p is density, k is thermal conductivity, and Cs is 
specific heat capacity of workpiece material in solid state. 

Boundary and initial conditions 


Energy transferred to the workpiece as heat input serves as thermal boundary 
condition on the top surface SI (Fig. 2.3). The bottom surface of the workpiece is 
assumed to be sufficiently away from the top surface so as to remain at its initial 
temperature throughout the grinding pass. The heat loss to the coolant on the surfaces S 1, 
S2, and S3 (Fig. 2.3) is modeled using convective boundary condition. Thus, the 
boundary conditions are: 


q=-kWT.n = -q,+hJ 

on 

S\,S2,S3 




where 






+ K ^0 

on 

S\ 

when 

t>0 

(2.3) 

q. = K ^0 

on 

SI, S3 




II 

on 

SA 

J 





Here, he is convective heat transfer coefficient. To is ambient temperature (i.e. dielectric 
temperature) and qwg is heat flux supplied to the workpiece due to grinding. 

The initial temperature T, can be taken as normal room temperature To of the dielectric 

in which workpiece is completely dipped. Thus, the initial condition becomes 

7; =r„ at t=0 (2-4) 

Galerkin Formulation 

The Galerkin weighted residual method [73] is used to convert the governing 
differential equation (Eq. 2.2) into algebraic equations by using appropriate weighting 
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functions. When the residue is made orthogonal to the weight function, the following 
integral equation is obtained. 




ar 

dt ' 


'd,,dT^ d ,,dT 
dx ox dz oz 


>W dxdz =0 


(2.5) 


where, W is weighting function and H is the domain of interest. In vector form, this can 
be written as 

( 2 . 6 ) 


f pC.^Wdxck -\ \y.{kVT)]Wclxct =0 
^ dt in'- ^ ^ 


which can be simplified to 

^^pC.^Wdxdz - |^[V. (iTArVr) ] dxdz + dxdz = 0 (2.7) 

The expression (2.7) is further simplified by using the divergence theorem, which gives 
rise to certain boundary terms. 

\^pC,^Wdxdz - ^^WkVT.n dT+^^VW.V{kT)dxdz =0 (2.8) 


where, n is outward normal to the boundary. 

Note that W=0 on S4 (the boundary on which T is specified). Using this condition as well 
as the boundary conditions of equation (2.3), the above expression becomes 

Ii + 12 + I3 =0 (2.9) 

Where 

A = f pC, — Wdxdz 
" dt 

I,=\vWN{kT)dxd. 

I-,=-lWq,dT + lh^WTdT 


Here, the boundary T is given by r=SlioS2<^S3 


( 2 . 11 ) 
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Finite Element Approximation 

In contrast to the traditional Galerkin weighted residual method where a global 
approximation is chosen for the field variable, in Galerkin finite element method the 
domain of interest is discretized into small elements [Fig. 2.4]. Then, the weighted 
residual integrals are evaluated over each element and finally the elemental expressions 
are assembled appropriately to get the global expressions. 





Fig 2.4 Finite element mesh of workpiece in cut-off grinding 

Different types of elements (triangles, rectangles, quadrilaterals, and curved 
triangles) are available for discretising the domain and choosing a suitable approximation 
for the primary variable. While selecting the element type and the corresponding shape 
functions, it must be ensured that the shape functions satisfy the following two mesh 
convergence criteria: 
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(1) completeness of the approximation in the limit of the element size approaching to 
zero, and 

(2) compatibility condition which states that the continuity across the interelement 
boundaries of the primary variable and its derivatives upto one order less than the 
highest order appearing in the weighted residual integral. 

There are quite a few 2-D elements which satisfy the above mesh convergence 
criteria. Amongst them, 8-noded quadrilateral isoparametric element seems to give better 
results. Therefore in the present work, 8-noded quadrilateral isoparametric element is 
used (Fig. 2.4) [13]. 

To obtain the finite element equation from expression (2.10), the domain is 
discretized into nem number of elements [Fig. 2.4], Further, over a typical element, the 
temperature is approximated using the unknown nodal values T;* and the known 
interpolation (or shape) functions Nj*. 

Thus, 

g 

T{x,z,t) = ^N:(x,z)T;{t) (2.12) 


where 

Nl jv;j (2.13 a) 

and 

{t-YAt; T- J,-] (2.13 b) 


In Galerkin’s method, the weighting functions are also expressed using the same shape 
functions. Therefore 

W(x,z) = 

t=l 

{vr}= & ==[£•] {r} 


(2.15) 
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where. 
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Similarly, 

^w}=[b'‘]{w^] 

Further, 


aT 

dt 




(2.16) 


(2.17) 


In the present work, the geometry of the workpiece is also approximated by the 
same shape functions: 

x={yV'}^ {.r"} , r = {z'} (2.18) 

The vectors (.V'} and [z‘^ contain the nodal values of x and z coordinates 
respectively. 

To express the boundary terms in terms of the nodal quantities, approximation 
over the boundary is needed. Note that, as the domain is divided into nem number of area 
elements, the boundary automatically gets divided into nbe number of line elements. The 
approximation for T, W, and qw over a typical boundary element can be written as: 

' .61 






T? 


(2.19) 


W = {n^} [ {iV* Nl iVj } 


w; 




( 2 . 20 ) 


?. ={A^‘f {?.*}={< V?} 




( 2 . 21 ) 
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Note that the vector {N**} contains 1-D quadratic shape functions, which are consistent 
with the 2-D shape functions Nj®. 

Finite Element Equations 

Substituting expressions (2.15-2.17) and (2.19-2.21), the integrals Ii, I 2 , 13 can be written 
as: 

pC,[n-]{n-} dxjz'fir 

e^\ 




{A'*HA'‘r‘/r.]{r}‘ (2,24) 

Here, Q® is the domain of typical area element while Tb is the domain of typical 
boundary element (Fig. 2.4). 

Then, the expression (2.9) becomes 

IW M{i"}+Ekr[0A-l {r}-X(»'rfG/*}+ 2 {rf(G*‘J{ 7-)‘=0 (2.25) 

where, 

[gc-] = \pc,\n-}{n-Y dxdi 
[g^:-] = \k[B-Y[B’]dxdz 

]pK‘\ = rfr I (2.26) 

{g/‘} = 

r* 

Here, [GC®] is elemental capacitance matrix, [GK‘]is element conductivity matrix, 
[G^^]is boundary conductivity matrix, and {g f is boundary flux vector. 

For the purpose of numencal evaluation of the integrals, the area integrals (over 
n®) are transformed to the natural coordinate (^, q) using the following transformation: 

I 1 

" J \(:'-)\J\d^dT} f2.27) 
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where, J is the determinant of the Jacobian matrix of transformation * 


[^1 = 


(2.28) 


dx dx 

dt] 
dz dz 

.aj 

The matrix [j] is evaluated using the transformation Eq. (2.18). The integrands of the 
integrals (2.26) also need to be expressed in terms of the natural coordinates ^ and tj. The 
choice of the natural coordinate system is dictated by the Gauss quadrature rule [73] used 
in the numerical evaluation of integrals over the element. 

Similarly the boundary integrals (over Tb) are transformed to the natural 
coordinate q using the transformation; 

= J(- )\J^\dg (2.29) 

r -1 

where, \Jf,\ is the Jacobian for the boundary element and is given by 


\Jh 




+ 


^ dx^ 


(2.30) 


Gauss-Legendre numerical integration scheme is used in evaluating the area as well as 
boundary integrals (2.26). The formula for area integral is 

I ^ (2.31) 

-/ -I 1=1 J=i 

while, that for the boundary integral is 

(2-32) 

-I '=1 

Here, n^, n^ and n,; denote the number of Gauss points in r\ and g directions 
respectively. Further, w; and wj denote the weights, and ^ i , r) j and q i denote the Gauss 
point coordinates. 

All the elemental matrices are evaluated using 3X3 Gauss quadrature for 
reasonably accurate evaluation of the integrals. Similarly, the elemental boundary 
integrals are evaluated using 3 point Gauss quadrature [73]. 



38 


Thermal Analysis of EDDG 


The final global element equations are obtained by assembling the elemental 
equations (2.26) at appropriate locations. The assembled finite element equation can be 
written as: 

[GGC] \t^ + [GGK] ^}^^GF} (2.33) 

Where [GGC] is global capacitance matrix due to grinding, [GGK] is global 
conductivity matrix due to grinding, [GGF] is global heat flux due to grinding, {T} is 

global temperature vector and \t | is time derivative of {T}. 

Equations (2.33) represent a set of ordinary differential equations in the variable T 
as function of time t. These equations are solved after application of implicit Finite 
Difference Method (FDM) discussed in subsequent section. 

Finite Difference Scheme 

The 9 family finite difference approximation for the time derivative ItI is used to 


convert the set of differential equations (2.33) into a set of algebraic equations: 

0<9<1 (2.34) 


where, {t}, 


and 


T are the temperature at step (i) and 

J 1+1 


step (i+I), 


and ( T 


are the derivatives of temperature with time at step ( i ) and step 

I+l 

(i+1) , and is the time interval of step (i+l). 

The resulting algebraic equations are: 

M,+i{^},+i = f = 0,l,2,3 n (2.35) 

where, 

M.+i =[C]+0 A/,.,1 [a:] (2.36) 

Wi,i+i = ^,+1 (| 9 {f},.,! + (i- 6 '){f},)+([c]-(i-^)a/,+, [ic:]){r}, (2.37) 

In the present case, the value of 9 is chosen to be 2/3 (Galerkin scheme). In each 
time interval At, the nodal temperatures at the new time level are calculated by using the 
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set of simultaneous algebraic equations (2.35). The solution thus marches in time, in 
steps of At until the desired final time is reached. In grinding. At is selected as time 
required to move down exactly by one element length of the grinding domain. 

2.2.3 Results and Discussion 

The temperature distribution on the top surface and along the depth are 
determined using FEM model and compared with the previously obtained experimental 
results [71] for the same machining parameters and material properties (table 2.1 and 
2.2). However, the feed value has been changed to 3.2 mm/s fi'om 3.38 mm/s for some 
computational convenience. Rectangular heat flux is assumed to be distributed over all 
elements of the top surface and it moves downward exactly by one element after every 
time step. Convergence study has revealed that the mesh of (19) X (22) isoparametric 
eight noded quadrilateral elements with 1337 nodes are adequate for this problem. 

Figures 2.5 show the temperature variation along depth as well as top surface 
temperature variation in dry cut-off operation at different grinding times. Figure 2.5a 
shows an agreement between the experimental and computed results obtained for 
different grinding times. The discrepancy between the computed and experimental results 
could be assigned to the difference in the feed value used (3.2 mm/s in place of 3.38 
mm/s used in reference [71]) and unknown grinding time for which experimental results 


1477 K[71] 


(a) 




(b) 


Fig. 2.5 Temperature variations with time during dry cut-off operation (a) depthvnse, (b) 
at the top surface, for U. =13.58 GPa and IC* =0.31 at Vf- 3.2 mm/s 





40 


Therma l Analysis of EDDG 


are reported. The workpiece top surface temperature was extrapolated using measured 
temperature at different depths below the top surface. The extrapolated top surface 
temperature is 1477 K for V,=3.38 mm/s. Using the present model, top surface 
temperature in similar conditions is determined as 1 124 K for Vf — 3.2 mm/s (Fig. 2.5b). 

The results for depthwise temperature variation are also compared at higher feed 
of 9.35 mm/s (Fig. 2.6). Here, results (Fig. 2.6a) are in close agreement if specific 
grinding energy Uo of work material and energy partition Rwg are considered feed 
dependent [71]. The values of Uo and Rwg for various feed rates are given in table 2.2. 
Theoretically calculated temperature values are higher as compared to the experimental 
values [71] if Uo and Rwg are considered feed independent (Fig. 2.6b). The top surface 
temperature in similar conditions is 1440 K (Fig. 2.6a where depth = 0) against the 
extrapolated value of 1144 K [71] if Uo and Rwg are considered feed independent. By 
choosing feed dependent Uo and Rwg (table 2.2), top surface temperature is obtained as 
855 K (Fig. 2.6b). One of the reasons for these discrepancies in the values of top surface 
temperature may be unknown time for which actual experiments were conducted. 


1600-1 



-(t=1.50s) 

-(t=1.71 s) 

• (t=1.93 s) 

■ -(1=2.1 4 s) 

- (t=2.35 s) 

~ (Experimental[71]) 



200 ^ 
0.000 


O.OOS 0.010 0.015 

Depth from top surhice(m) 


0.020 



(t?«2,14») 
- (tw24S#) 


0.000 0J)0S 0J>10 04)1S 0J20 

D«ptti from top surf«co(m) 


(a) 


(b) 


Fig 2.6 Depthwise temperature Variations with time during dry cut-off operation at V, »9.35 
(a) Rwg =0.27 and U, =6.2 GPa (a) R,* =0J1 and U, =13.58 GPa 
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From above comparison it is evident that accurate prediction of temperature can 
be made if accurate values of Uo and Rwg are known in advance. 

Figures 2.5 and 2.6 show that as one moves down from the top surface the 
temperature decreases, reaches a minimum value, and then increases. This phenomenon 
is observed because of higher thermal front speed than the rate of down feed into the 
workpiece and consideration of bottom surface always at room temperature. 

Temperature distribution 

Effects of different feeds, energy partitions, specific energy, different dielectrics 
(as coolant) and diamond grinding wheel on temperature distribution are discussed 
below. 

Feed 

Figure. 2.7 shows the calculated top surface temperature distribution obtained at 
different feeds. Here, it is observed that the temperatures are higher at higher feeds if Uo 
and Rwg are considered feed independent (Fig. 2.7a). In similar conditions top surface 
temperatures are lower (Fig. 2.7b) at higher feeds if Uo and Rwg are considered feed 
dependent (table 2.2) because Uo and Rwg both decrease with increase in Vf for the 
constant heat flux (Eq. 2. 1). 

Fig 2.8 shows the depthwise change in temperature. It is also observed that when 
Rwg and Uo are taken feed independent then the temperature gradient are higher for higher 
feed (Fig. 2.8a). But, trend becomes reverse if feed dependent values of Rwg and Uo are 
considered (Fig. 2.8b). 

Further, accurate prediction of temperature is possible only when sufficient data 
are available for variation of Uo and Rwg with different feeds for different grinding wheel 
and workpiece combination. At higher feed, initially a high temperature is generated due 
to higher values of Uo, Rwg, and Vf (Eq. 2.1) but after some grinding time Uo and Rwg 
start to decrease due to temperature softening effect 
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Fig 2.7 Top surface temperature variations during dry cut-off operation using Alj Oj 
grinding wheel for different feed rates after 2.5 s of grinding (a) =* 13.58 GPa and 

Rw ,5 = 0.31, (b) Uo and are feed dependent and selected from table 2.2. 



(a) W 


Fig 2.8 Depthwise temperature variations during dry cut-off operation using Alj O 3 
grinding wheel for different Vf after 2.5 s of grinding (a) U, = 13.58 and R., =» 0.31, 
(b) U, and R„g are feed dependent and selected from table 2 . 2 . 




Thermal Analysis of EDDG 


43 


Energy Partition 

Figure. 2.9 shows the temperature distribution obtained at different energy 
partition values. Here, it is observed that the temperatures are higher at higher Rwg 
because of higher heat flux supplied to the workpiece which is calculated by using Eq. 
2 . 1 . 



0.000 0.006 0.010 0.018 0.020 
Depth from top surface (m) 



0.000 0.009 0.010 0.015 0.020 

Distance along x<axie (m) 


Fig 2.9 Temperature variations during dry cut-off grinding using alumina wheel for 
different R*, at U, = 13.58 GPa and Vf= 6.4 mm/s after 2.5 s of grinding (a) along depth 
(b) top .surface. 


Specific Energy 

Figure. 2.10 shows the temperature distribution obtained at different specific 
energy values. Here, it is observed that temperatures are higher at higher Uo because of 
higher heat flux supplied to the workpiece, which is calculated by using Eq. 2.1. 
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Fig 2.10 Temperature variations during dry cut-off grinding using alumina wheel for 
different U. at = 0.27 and Vf= 6.4 mm/s after i5 s of grinding (a) along depth, (b) 
top surface. 
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Dielectrics 

Figure 2.11 shows the efifect of different dielectrics on top surface temperature 
variations with alumina wheel at Vf =4.0 mm/s after 2.5 s of grinding. The data related to 
the workpiece are taken fi-om table 2.1. The top surface temperature drops to half of its 
value when a coolant having cooling capacity of 25000 W/m K is used in place of dry 
grinding. Curvature in top surface temperature plot is observed at two ends due to 
convection at surfaces S2 and S3 (Fig. 2.3). 
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Fig 2.11 Variation of top surface temperature during cut-off grinding using alumina 
wheel for different h, at Vf =4.0 mm/s, U„ =13.58 GPa, and R*, =0.31 after 2.5s of 
grinding 

Diamond wheel 

Figure 2.12a depicts the variation of top surface temperature with time during 
diamond cut-off grinding ofHSS workpiece. Coolant is also applied during the process. 
The material property data are given in table 2.3. The temperature is found to be lower 
with diamond wheel in comparison to alumina wheel (Fig. 2. 12b) because of lower value 
ofRwg(0.08) used for the case of diamond cut-off grinding. Further thermal diffusivity of 
diamond is much higher (due to higher thermal conductivity) and hence major portion of 
heat generated during grinding is taken away by the grinding wheel. No information is 
available in literature for actual value of Rwg during diamond cut-off grinding. The value 
of Rwg varies fi'om 0.08 to 0.18 with different feed rates if special graded alumina 
abrasives are used for cut-off grinding of AISI 1020 steel workpiece [71], The value of 
Rwg for grinding of HSS with diamond wheel has been assumed to be approximately 
equal to the value of Rwg for grinding of AISI 1 020 steel with special graded alumina 
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wheel. Wheel clogging is neglected but in actual conditions this is a major problem in 
diamond grinding of steel. 



Fig 2.12 Variation of top surface temperature with time during cut-off grinding 
using (a) diamond wheel, (b) alumina wheel at Vf =10.29 mm/s, U„ =10.55 GPa, 
and ht =10000.0 W/m^K (kerosene) 


2.3 Thermal Analysis of EDM 

EDM is another constituent process of EDDG, which is basically a thermal process. 
Large literature is available regarding thermo- mathematical modeling of EDM process, 
which is summarized in section 1.5.3. Here, the following aspects of EDM are discussed 
first; spark radius at workpiece-flux contact surface, quantitative value of flux, energy 
partition of flux, and latent heat consideration for phase change. 

2.3.1 Mathematical Modeling 

Spark radius (R) 

It is found that electrode material influences the discharge channel diameter during EDM. 
However, quantification/measurement of spark radius are extremely difficult due to very 
high pulse frequencies. Erden [74] suggested an empirical relationship to evaluate radius 
of a spark for the selected pairs of electrodes and dielectric. For a rectangular pulse, the 
spark radius R(t) as a function of time (t) is written as: 

Rit) = KQ’”t’' (2.38) 

where Q is discharge power, and trt, n and K are empirical constants. These constants are 
defined in terms of experimental coefficients L, M and N determined experimentally for 
different electrode materials; 
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_■ ^ m = QSN n = N (/ is discharge length) (2.39) 

Snoeys and Van Dijck [55] varied heat source diameter and studied its effects on various 
responses. Shankar et al. [70] have found noncylindrical spark shape with different radii 
at different cross sections across the spark length based on their integrated model wherein 
temperature distributions in the electrodes and dielectric are analyzed together. Pandey 
and Jilani [62] have proposed the following model for calculation of discharge radius R: 


^ EaR -1 
r, =-^tan 


k n 


05 


Aat 


-,05 


(2.40) 


where, Tb is boiling temperature, Eo is energy density, a is thermal diffiisivity, t is the 
time at which spark radius is calculated, and k is thermal conductivity of workpiece 
material. But, the application of this expression is limited, as it is valid only for a few 
cases of selected electrode pairs and dielectric. Literature survey indicates that although 
importance of spark radius during EDM process has been acknowledged, there is a 
scarcity of comprehensive work on the evaluation of its shape and size. In the present 
work the radius of a spark is calculated based on the results of reference [70], 

Heat flux due to a single spark 

Most researchers [63,74,75] have considered uniform disc heat source distribution 
within a spark. This assumption is far from reality. This fact is evidenced from the actual 
shape of a crater formed during EDM. In the present work, a Gaussian heat flux 
distribution (Appendix A) is assumed. The heat flux qw* (r) at a radius (r) can be 
calculated if maximum intensity go at the axis of the spark and its radius (R) are known. 


= - 4.5 



W 


(2.41) 


If it is assumed that total energy of each pulse is to be used only by one spark, then the 
above equation can be written as follows (Appendix -A); 




4.45 i?,, U, I 


'CXp 




(2.42) 


where Ub is breakdown (discharge) voltage, I is current and Rws is the energy 
paitition(fraction of input heat going into the workpiece) due to EDM. 
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Energy partition (R^ in EDM 

Another important parameter required for computational analysis of EDM process 
is heat distribution between cathode, anode, and dielectric (neglecting radiation). The 
material properties of individual electrodes determine the fraction of input heat going into 
the electrodes. Theoretically, steel anode gets lower share of the input power than copper 
due to its lower thermal diffrisivity. 

Researchers [65-67] assumed that different fraction of the total discharge power 
(= Ub I) is lost to cathode (18%), anode (8%) and dielectric (remaining 74 %). The 
fraction that goes to cathode and anode is independent of current and pulse-time (tp). The 
validity of this assumption was verified by DiBitonto et al. [65] by showing the 
agreement between theoretical and experimental optimum pulse times with same erosion 
rate. Shankar et al. [70] have theoretically calculated that about 40-45 % of the heat input 
is absorbed by the workpiece made of 5 Cr die steel. They have calculated this value 
using water as dielectric and copper as tool electrode. But nobody has given a 
comprehensive way to calculate theoretically or determine experimentally the value of 
Rw, during EDM process. In the present work, R^ is taken as 0.08. 

Latent heat for phase change 

The latent heat required for phase change is accounted for by modifying the 
expression for specific heat of workpiece (Cs). Considering the enthalpy before and after 
phase change, latent heat of melting (Lm) and latent heat of evaporation (Lv) are 
incorporated into modified specific heats Ci, and Cv in the vicinity of melting temperature 
Tm and boiling temperature Tv, respectively [69]: 



forT„-^T<T<T„+^T 

(2.43) 


forT^-^<T<T, + ^T 

(2.44) 


where Ci and Cv are specific heats of workpiece material m liquid and vapor states, 
respectively. 
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2.3.2 Finite Element Formulation 

Governing equation 

Heating of workpiece due to a single spark (Gaussian distributed) is assumed to be 
axisymmetric i e — — 0 Therefore, temperature field is goyemed by thermal diffusion 

' ' ' de 


differential equation 


PCs 


dt 


I d ,, dT . d ,, dT . 
r dr dr dz dz 


(2.45) 


where, r and z are coordinate axes (Fig. 2.13). 

Boundary and initial conditions 

A small cylindrical portion of the workpiece around a spark is used as the domain. 
Energy transferred to the workpiece as heat input serves as thermal boundary condition 
on the top surface Ti (Fig. 2.13). Heat loss to the coolant on the surface Fi is modeled 
using the convective boundary condition. Further, the boundaries r 2 and Fs are at such a 
large distance that there is no heat transfer across F 2 and F 3 (Fig. 2.13). Similarly, there 
is no heat transfer across the axis of axisymmetry i.e. across F 4 (Fig. 2.13). 

The heat transfer across any boundary is given by 


q = -V T.n = + hT 


where, 

for r<R 
=KTo for r>R 
= 0 for off - time \ 
h =0 for r < R 
= -4 for r>R 
= 0 for off - time 


for t > 0 on Fi 


(2.46) 


P. on F 2 , Tj and r 4 (r = 0) (2.47) 

Here, qws is quantity of heat flux entering into the workpiece, R is radius of spark, Rw* 
is discharge energy partition fraction and n is normal to the surface. 
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The initial temperature T, can be taken as temperature of the dielectric in which 
workpiece is completely dipped. 

^ ~ ^ at t=0 (2.48) 



Fig 2,13 Thermal model of EDM 

Galerkin formulation 

Here also, Galerkin weighted residual method [73] is used to convert the 
governing differential equation into algebraic equations by using approximate weighting 
functions. When the residue is made orthogonal to weight functions, the following 
integral equation is obtained; 

f pC,—W27trdrdz 

^ dt ^r2\jdr dr dz\ dz )_ 

The expression (2.49) is simplified using divergence theorem, which gives rise to certain 
boundary terms; 

\MV?^rdnk + f \ 

^ dt ki^dr dr dz dzj 

Using the boundary conditions (2.46), the above expression (2.50) becomes 
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Finite element approximation 

To obtain the finite element equation from expression (2.52), the domain is 
discretized into nent number of elements (Fig. 2.14). Further, over a typical element, the 
temperature is approximated using the unknown nodal values T* and the known 
interpolation (or shape) functions Ni®. In the present work, 8-noded quadrilateral 
isoparametric element is used (Fig. 2.14). 



Fig 2,14 Finite element mesh in EDM domain 
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Thus, over a typical element 

g 

= Y^N:(r.z)T:{t) ={,V'}"{r} (:.J3) 

|S*l 

where 

{vV'}^={,v,' n: n; ;v; } (2.54) 

For Galerkin's method; 

P 

^V(r,z) = X‘y:(r,:)^v; = {y^')"{lf'}={(r‘}"{iV'} (2.55) 

t = t 

Then, the expression for derivatives of W and T becomes: 
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(2.56) 
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= k]{r} 

2L^{n-y ■] 


(2.57) 


(2.58) 


(2.59) 

(2.60) 


Here, the nodal values of W are known but arbitrary. The geometry is also approximated 
by the same shape functions: 

r^N’Yir-} --{iV-l'h-) < 2 ^ 61 ) 

The vectors {r*} and {r'} contain the nodal values of r and z coordinates 
respectively. 

^ r * ’ IT T’WN 

iBTtcfhr s't . .i V. ff-A 
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To express the boundary conditions in terms of the nodal quantities, 
approximation over the boundary is needed. Note that, as the domain is divided into nem 
number of elements, the boundary automatically gets divided into nbe number of 
elements. Now approximation for W, T, and qw over a typical boundary element can be 
written as: 


w = {n‘’Y 


K ^ 3 } 








(2.62) 


(2.63) 





(2.64) 


Note that the vector (N*’} contains 1-D quadratic shape functions, which are consistent 
with the 2-D shape functions Nj*. 


Finite element equations 


Substituting expressions (2.56-2.60) and (2.62-2.64), the integrals hM and I, can be 
written as 


A =|W'i pC,{N-lN-Ydrdz)l^TA 

nem / 

e=l / \ J 


nbe . /■ 

WN‘Nr. U:]] MN‘^r.){rr (2.67) 

Here, Q* is the domain of typical area element while T, is the domain of typical boundary 
element. Then the expression (2.51) becomes 



( 2 . 68 ) 
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where. 


[5C']= jpC, rdrdz 

[5A:']= \k[B^Y[B‘\rdrd 2 


= rdV 


(2.69) 


Here, is elemental capacitance matrix, [SST*] is element conductivity matrix, 

is boundary element conductivity matrix, and (iSy* is boundary element flux 


vector. 

When the elemental quantities of Eqs. (2.69) are assembled, the following global 
equations are obtained. 

[aye] |r| + [G 5 '^] ^}={ gsf } (2.70) 

Where [GSC] is global capacitance matrix due to spark, [GSK] is global conductivity 
matrix due to spark, [GSF] is global heat flux due to spark, {T} is global temperature 

f‘l 

vector, and derivative of {T}. 

Equations (2.70) represent a set of ordinary differential equations in the variable T 
as a function of time t. These equations are solved after application of implicit Finite 
Difference Method (FDM). 

Finite difference scheme 

The finite difference scheme used is the same as discussed in section 2.2.2 with the 
differences that [C], [K] and (F) are replaced by [GSC], [GSK] and {GSF} respectively. 
Here, two-time steps, Ati for the on-time and Ati for the off-time of EDM pulse power 
supply, are selected. 
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2.3.3 Results and Discussion 


In order to assess the accuracy and sensitiveness of the EDM model, a problem 
involving variable heat flux profile over an axisymmetric zone is considered. Results are 
compared with the known theoretical results of Shankar et al. [70], Since the heat flux 
profile is not of a standard configuration in Shankar’s model, the results are obtained 
using two commonly used shapes: Circular disc heat source and Gaussian distributed heat 
source. The effect of heat flux is assumed [70] to be significant only upto a distance of 
4R from spark center along radial direction as well as depthwise direction. The material 
properties are given in table 2.4 and process parameter in table 2.5. The heat flux is 
distributed over four elements along radial direction at workpiece top surface. From 
convergence point of view, it has been found that 256 elements with 833 nodes are 
adequate for the specified accuracy. Top surface temperature variation after 60 |4S (after 
one spark) along radial direction is shown in Fig. 2.15. The results of present model with 
circular disc heat source are in closer agreement with the results of Shankar et al. [70] 
than the one with Gaussian heat source. But temperature distribution trend with Gaussian 
heat flux is closer to the shape of the crater formed at the workpiece surface by a spark. 


—Shankar et al (70) 

• - - Present model-Gaussian heat source 
Present modei-Oisc heat source 



40 60 80 100 120 140 160 180 200 

Radial Distance from centre ( M' m ) 


Fig 2.15 Comparison of variation of top surface temperature in 5 Cr die steel 
workpiece during EDM with heat flux of different profdes and the results of 

Shankar et al [70]. t=60 ps, t,. =€0 ps, t, =120 ps, without considering phase 
change and on-/off-time effect. 
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Effects of On-/ofF-time, duty cycle, energy partition and current on temperature 
distribution are discussed below. 

On-/off-time 

Figure 2.16 shows the top surface temperature variation considering the 
effects of on-/off-time and phase change (change of state from solid to liquid and liquid 
to vapour). The material properties are given in table 2.4 and process parameters are 
given in table 2.5. Further, to study the effect of on-off time on temperature distribution 
the time selected is 60 p.s which falls in the off period of the pulse. The temperature is 
calculated using equation 2.70 and taking the time step At = 30 ps (= on-time). 
Computational results can be obtained after a time interval, which is an integral multiple 
of on-time. Hence, the results are given after 60 ps. 



—without phase change and on-off both 
» --with phase change only 
—with on-off without phase change 
- with phase change and on-off both 


} 25 50 75 100 125 150 

Radial Distance from centre ( jim ) 


Fig 2,16 Top surface temperature variation in 5 Cr die steel workpiece during EDM computed 
with Gaussian heat flux profile including the effect of phase change and on-/off- time both 
after 60 ps for t« =30 ps, tp =120 ps 


Obviously, higher peak temperature is observed if the effects of off-time of a 
pulse and phase change in the work material are not considered. This is so because it is 
assumed that the heat is supplied over the entire pulse period (off-time = 0), and further 
no heat is utilized for the change of phase from solid to liquid and then liquid to vapour. 
Further, reduction in temperature is observed if only phase change effect is considered 
because of the fact that some (latent) heat is utilized for phase change at constant 
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temperature. Therefore, the temperature rise is less in comparison to the previous case. If 
both the effects (on-/off-time and phase change) are taken into consideration, the peak 
temperature is still lower. This is because of the fact that heat is supplied only during the 
on-time of the pulse, and no heat is supplied during the off-time. Also, some fraction of 
the supplied heat is also spent on phase change. Hence, the overall temperature rise is 
much lower. This case is closer to real life situation. 

Duty cycle 

Figure 2.17 shows the effect of duty cycle (which is defined as the ratio of on- 
time to pulse time) on top surface temperature distribution in HSS workpiece after one 
pulse time. The duty cycle is changed by changing the on-time of the power supply by 
keeping the off-time same. The material properties are given in table 2.3 and process 
parameters are given in table 2.6. The curves show that higher temperature is obtained 
with higher value of duty cycle with the same off-time. About 10 % increase in peak 
temperature is observed by increasing the cycle time by 25%. 

From above, it can be concluded that the present model simulates EDM process 
satisfactorily. It accounts many realistic facts viz., consideration of heat loss to the 
dielectric medium, on-/off-time consideration, energy partition to workpiece, duty cycle 
current, and consideration of phase change (Fig. 2.15-Fig. 2.20). 


duty cycle>=40 % 
- duty cycle«60 % 
duty cyclesSO % 



25 50 75 100 125 160 175 200 
Distance along radius ( pjn ) 


(a) 


I* - • • duty cycle«40 % 
duty cycl««50 % 
- duty cycle»60 % 



25 50 75 100 125 ISO 175 200 
Depth from top surface ( pm ) 


(b) 


Fig 2.17 Variation of temperature with duty cycle (a) at top surface, (b) along depth, in 
HSS workpiece during EDM computed with Gaussian heat flux after one pulse time for 

toflplOO ^is. 





Depth from top surface (jun ) 


Fig 2.18 Depthwise temperature variations in HSS workpiece during EDM 
with Gaussian heat flux for t« =100 ps, t ,=200 ps after 100 ps at radial 
distance of (a) r = 0, (b) r = R, (c) r = 2R 








58 


Thermal Analysis of EDDG 


Energy partition (R ws) 

In the present FEA model prediction of temperature is governed by another 
important parameter, i.e. energy partition (Rws). The thermal properties of workpiece 
electrode are given in table 2.3. Theoretically, electrode having lower value of thermal 
diffusivity (k/pC,) receives lower share of heat energy than the one having higher value. 
Fig. 2.18 depicts the effect of Rw. on depthwise temperature distribution at different 
radial distances from the center of spark (r = 0, r = R and r = 2R) after 100 (as. Process 
parameters are given in table 2.6. Higher temperature is observed at higher R^,,, which is 
expected because it gives more amount of heat entering into the workpiece. Trends are 
similar at all radial distances except at r =2R which is lying outside spark radius. The 
curves show that the temperature gradients are high along the depth of the workpiece if 
radial distance is small. 

Figure 2.19 shows the top surface temperature variation along the radius for 
different values of Rw, after 100 (as. The top surface temperature at the center of the spark 
is also three times high when Rws is increased by four fold. The temperature gradients at 
the top surface (within the spark radius) are high but comparatively less than the 
depthwise temperature gradient (Fig. 2. 18a). 



Fig 2.19 Top surface temperature variation in HSS workpiece during EDM 
with Gaussian heat Hux for different values of partition ratio (R,^ at t«. -100 
ps, tt,=200 ps after t=100 ps. 
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Current 

Figure 2.20 shows the effect of current on temperature distribution at top surface as 
well as along the depth after 100 jis. The workpiece (HSS) material properties are given 
in table 2.3 and the process parameters are given in table 2.6. It is observed that increase 
in current leads to increase in peak temperature due to increase in heat flux (Eq. 2.44). 
The temperature gradient along depth is found to be steeper than temperature gradient 
along radial distance from center. 

It is also noted that the temperature gradient is very high upto 30 |j,m depth from top 
surface for all current values. Hence, from thermal stresses point of view, this thin 
surface layer is most critical. After this depth, temperature gradient gradually decreases. 
Almost no change in temperature at top surface as well as along depth is observed 
beyond the distance of one spark radius. 




Depth from top eurfece^m) 


Fig 2.20 Temperature variation in HSS workpiece during EDM with Gaussian heat 
flux for different discharge current at t,. =100 ps, duty factor=50% after t=l00 (a) 

along depth from top surface (b) along radial distance from centre 


2.4 Temperature distribution in EDDG process 
2.4.1 Theory of Superposition 

Since, temperature independent properties of workpiece material are considered, 
the heat conduction equation is linear, hence the sum of known solutions also satisfies the 
heat conduction equation. Under these circumstances the combined temperature variation 
with time can also be obtained by addition of known solutions of simpler component 
problem [76], According to this method different boundary temperature distributions with 
different initial conditions can be added to find final temperature distribution. Hence, the 
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resulting temperature (TT) due to grinding and EDM both is obtained by superposition of 
the two temperatures as follows; 

TT=To+GTR+EDMTR (2-71) 

Where, GTR is temperature rise due to grinding and EDMTR is temperature rise due to 

EDM. 

2.4.2 Results and Discussion 

The EDDG in cut-off grinding mode is used for a HSS workpiece (6 mm X 6 mm 
X 20 mm) (Fig. 2.21a and table 2.3) with a grinding wheel having its width (6.5 mm) 
greater than the width of the workpiece (6 mm). Top surface discretization and locations 
of sparks (points •) on the workpiece surface are shown in Fig. 2.21b. Superposition 
technique is used to obtain temperature distribution in the workpiece due to EDDG, 
which includes the effect of diamond cut-off grinding as well as EDM on the rise in 
temperature at a point (a node). 

First, 3-D temperature distribution in the workpiece due to cut-off grinding is 
obtained by using 2-D thermal analysis of cut-off grinding (Fig. 2.3-2.4). Let us consider 
a plane Y 1 Y 2 Y 3 Y 4 parallel to y-z plane (Fig. 2.21a). It is assumed that the temperature 
along the width of the workpiece (y-axis) at a specified depth in the given plane, say 
Y1Y2Y3Y4, is same during cut-off grinding. For example, temperature at any point along 
the line Y 1 -Y 2 on the top surface of the workpiece will be same during cut-off grinding. 
In the same way, temperature at any point along the line PPi (Fig. 2.21a) will also be 
same but of different magnitude as compared to that along Y1-Y2 in the plane Y1Y2Y3Y4. 
Using this procedure, 3-D nodal temperatures are obtained at any specified grinding time 
which is an integral multiple of time required for grinding wheel to move downward by 
one element length ( 0.5 mm). 

During EDM, sparking takes place at a point of least resistance [ 6 ], which 
depends on surface roughness of both the electrodes. Chances of occurrence of next spark 
at the location in the next cycle are comparatively less. Until spark does not reoccur at the 
same location, there will be no heat input at that point (or node). Keeping this 



Thermal Analysis of EDDG 



Fig. 2.21 (a) 3-D domain of workpiece (b) Descretization of workpiece top 
surface used in superposition 
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characteristic of EDM in mind, heat input to the workpiece by sparks has been simulated 
as follows. 

It has been assumed that spark strikes at the workpiece top surface at selected 
points (Fig. 2.21b). It is considered that only one spark takes place at a time (say, at point 
1) and next spark at the same location is repeated in sequence after sparking at all other 
selected locations (2 to 9 in Fig. 2.21b) are over. This pattern has been assumed to 
simplify the problem as well as to match the grinding time for one element movement in 
downward direction (in FEM computation) with total time of sparking in that duration. 
The grinding time for one element movement in downward direction is completely 
matched with total time of sparking at all other 9 locations. For this matching of time, 24 
times sparking is needed at each of the nine spark locations with on time of 75 jis and off 
time of 150 ps. In this way, at each spark location, heat will be supplied for 75 ps and 
again it will be repeated after 1950 ps (sum of the cycle time of 8 sparks of 8 locations 
plus off time of the spark at present location). 

After several numerical experiments it has been found that the axisymmetric 3-D 
domain influenced by the heat of a spark is upto 8R distance along radius as well as depth 
from the center of the spark after EDDG for 0.05 s (Fig. 2.22). This range of 8R may 
change with the time of machining (say, 0.05 s in this case). It is seen (Fig. 2.22) that the 
temperature achieved at 8R after 0.05 s is almost room temperature. Further, this whole 



Distance along radius ( *4 m ) 


1200 



3004 — — n — ^ 

0 200 400 «00 800 1000 

Deptti from top suitace ( pm ) 


Fig 2.22 Variation of temperature (a) at top surface along radius (b) along depth at r«0, 
of two different size of HSS workpiece domain after 0.05 s of machining time using EDM 
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axisymmetric domain is found to lie within 4-eIements length along x- and y-axis, and 
two-element length along z-axis of grinding domain. Nodal temperatures of the 
axisymmetric domain due to sparks, which lie on the nodes of grinding domain, are 
scalarlly added to the temperature rise due to grinding to get combined temperature (i.e. 
the change in nodal temperature due to EDDG process as a whole). 

The rectangular parallelepiped domain of 16R X 16R X 8R around a spark could 
be reasonably considered as a domain independent of the effect of heat from the 
neighboring sparks (assuming that the temperature distribution is evaluated after 0.05 s of 
grinding). Nine such solid blocks are assumed to be fitted within a single solid block of 
dimension 48R length along x-axis, 48R width along y-axis, and 8R depth from top 
surface opposite to z-axis. 

Due to spark the temperature is affected only in certain small domain around the 
spark. If the subsequent spark occurs outside the zone of influence of the first spark then 
it will not affect the temperature distribution given by the proposed model. Thus, even if 
spark locations are chosen randomly the temperature distribution will not change as long 
as subsequent spark locations are outside the zone of influence of the current spark. 
However, if the subsequent spark occurs within the zone of influence of current sparks, 
the temperature distribution will certainly gets affected but then the problem will become 
3-D. To avoid this complexity it is assumed that the subsequent spark occurs outside the 
zone of influence of current spark. 

Fig. 2.23(a) shows the top surface temperature distribution after superposition of 
grinding and EDM temperature distributions. Process parameters are given in table 2.6. 
The figure shows that the peak temperatures are observed at the nine spark locations. The 
contour plot of isotherms is shown in Fig. 2.23(b). It is seen under microscope that the 
EDM’d surfaces have craters all around it. These very high temperature peaks (Fig. 
2.23a) are responsible for the craters on the EDM’d surface. However, during EDDG, 
most of these craters would vanish during grinding part of the EDDG process. It is 
evidenced by the micrograph of EDDG’d surfaces given in reference [29]. 
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Fig 2^ <l«Wbution in HSS workpiece during EDDG process for t«,-75 fis, 

t,fl=150 US after t=O.05s (a) 3-D surface plots (b) Contour plots 
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Figure 24(a) shows the temperature distribution at workpiece top surface due to 
EDDG process in X-X plane (parallel to x-axis) at y=1.0 mm (Fig. 2.21b). As expected, a 
sharp increase in temperature is observed around spark locations. Top surface 
temperature of about 400 K is calculated due to grinding when no sparking is allowed to 



(a) 



(b) 

Fig 2.24 Top surface temperature distribution in HSS workpiece during EDDG after 
t^.OSs. (a) along x-axis at y=1.0 mm (along X-X in Fig. 2.21 b), (b) along y-axis at x=3.5 
mm (along Y-Y in Fig. 2.21 b). tttii=75 fis, tofr=150 ps 





66 


Thermal Analysis of EDDG 


take place by switching off the power supply. The peak temperature due to grinding and 
EDM both at the spark location is calculated about 1206 K after 0.05 s. Temperature 
gradient is found to be very high near the spark locations. Similarly top surface 
temperature variation along Y-Y (parallel to y-axis) at x = 3.5 mm (Fig. 2.2 1 b) is shown 
in Fig. 2.24(b). 

Depthwise variation of temperature at locations A, B, C and D (Fig. 2.21b) in the 
workpiece due to EDDG is shown in Fig. 2.25. It is noted that temperature gradient is 
very high within a thin surface layer of 0.63 mm. Beyond this depth temperature gradient 
variation is almost zero. This shows that temperature-softening effect is limited to a thin 
layer near the spark flux. 

Effects of duty cycle, on-time, current, energy partition, feed velocity, and time of 
machining are discussed below (Fig. 2.26-2.34). The process parameters other than the 
parameter under study are taken from table 2.6. 

Duty Cycle 

Figure 2.26a and 2.26b shows the effect of duty cycle on top surface temperature 
distribution in the workplace due to EDDG process. The process parameters are given in 
table 2.6. Duty cycle is changed by changing off-time and keeping on-time constant. No 
significant change in top surface temperature is noticed due to variation in duty cycle. 
This is because, after the on-time duration the heat supply at each of the spark location 
resumes after a long period (which is the sum of off-time duration of the present location 
plus the pulse time of other eight spark locations). Hence, the temperature of the EDM 
domain comes down to the same level for all duty cycles. 

Figure 2.26c shows the effect of duty cycle on depthwise temperature variation 
for the same set of conditions as used in Fig. 2.25. Very small change in temperature is 
observed upto 175 pm depth. This shows that the temperature-softening zone at the top 
surface does not change with the duty cycles. 

On time 

igure 2.27 shows the effect of on-time on top surface and depthwise temperature 
bution in the workpiece due to EDDG process. The process parameters are same as 
m table 2.6. Sigmficant change in top surface temperature is noticed due to variation in 
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(c) (d) 

Fig 2.25 Dcpthwise temperature distribution in HSS woriipicce during EDDG process after t=0.05s 
at (a) x=1.0 mm and y=1.0 mm (Point A), (b) x=1.5 mm and y=1.0 mm (Point B), (c) i=1.0 mm and 
y=0.5 mm (Point C), (d) r“1.5 mm and y=0.5 mm (Point D). (Fig 2.21b) 
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Fig. 2.26 Effect of duty cycle on temperature distribution (a) on top surface along x-axis at y-1.0 
mm (along X-X) (b) on top surface along y-aiis at r=3.5 mm (along Y-Y) (c) along depth at x»L0 
mm and y=1.0 mm (Point A) (Fig. 2.21b) in HSS workpiece during EDDG process after t*0.05s 
for Vf = 10.29 mm/s, t„.=150ps 
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Distance along x-axis (m) 


(a) 



(b) 

Fig 2,27 Effect of on time on temperature distribution (a) at top surface along x-axis at 
y=1.0 mm (along X-X) (b) along depth at x=L0 nun and y=L0 mm (Point A) (Fig. 2.21b) in 
HSS workpiece during EDDG process after 0.05 


on-time. This is so because larger amount of heat is supplied at each of the spark location 
during longer on-time than the shorter on-time pulses. 

Feed 

Figure 2.28 (a) shows the top surface temperature distribution in X-X plane at different 
feed rates. The process parameters are same as in table 2.6. At every spark location, the 
temperature at the center of a spark remains the same, as there is no change in electrical 
power input related to EDM. Secondly, after one spark has occurred, the heat supply by 
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the next spark at the same location resumes after a long period (which is the sum of ofF- 
time duration of the present location plus the pulse time of other eight spark locations). 
However, number of sparks at a spark location with increase in feed may decrease which 
may lead to higher contribution by grinding in EDDG. 



(a) 



(b) 

Fig 2,28 Effect of feed velocity on top surface temperature distribution in HSS 
workpiece during EDDG process after 0.05 s (a) along x-axis at y-1.0 mm (along X- 
X in Fig. 2.21 b) (b) along y-aiis at x=3.S mm (along Y-Y in Fig. 2.21 b). 


Away firom the spark center, the effect of feed is noticed. Higher temperature is 
observed for lower feed and vice versa (Fig. 2.28a). This is because, the contact duration 
of grinding wheel with workpiece is more at the lower feed. Similar trend is observed 
along y-axis at x=3.5 mm (Fig. 2.28(b)). Figure 2.29 shows the effect of feed on 
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depthwise temperature distribution. Increase in temperature is noticed due to decrease in 
feed. 



Depth from top surface (^m) 



Depth from top surface (m) 


(b) 


Fig 2.29 Effect of feed velocity on temperature distribution along depth in HSS workpiece 
during EDDG after 0.05 s at (a) j?=1.0 mm and y=1.0 mm (Point A Fig. 2.21 b) (b) x=1.0 
mm and y=0.5 mm (Point C Fig. 2.21 b). 


Current 

Effect of current on top surface temperature distribution is shown in Fig. 2.30(a) 
and 2.30(b). A significant change in the peak temperature is observed at spark locations. 
This is due to the fact that heat flux increases with the increase in current (Eq. 2.44). 
About five fold increase in temperature is observed when the current is increases to three 
times. The peak temperature is about 5900 K when the current is 24 A along X-X at 
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y=1.0 mm (Fig. 2.21b). Fig. 2.30c shows the effect of current on depthwise tem.perature 
distribution at point A of Fig. 2.21b. The effect is observed upto 500 ^im. 






0.000 0.001 0.002 0.003 0.004 
Depth from top surface(m) 


0.005 


Sonf distribution (a) top surface along x-axis at y-1.0 mm 
mm and y=1.0 mm (Pomt A Fig. Z21 b) in HSS woricpiece during EDDG after 0.05 s 
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Energy Partition due to EDM (R^ 

Effect of energy partition on top surface temperature distribution is shown in Figs. 
2.31(a) and 2.31(b). A significant change in the peak temperature is observed at spark 
locations. This is due to the fact that heat flux increases with R,v* (Eq. 2.44). The peak 
temperature is about 7200 K when the Rw, is 0.40 along X-X at y=1.0 mm (Fig. 2.21b). 
Fig. 2.32 shows the effect of Rws on depthwise temperature distribution at point A of Fig. 
2.21b. A sharp decrease in temperature gradient is observed upto 250 pm depth and then 
decreases slowly upto 500 pm before approaching to grinding temperature. This shows 
that the critical zone for high thermal stresses (because of high temperature gradient) lies 
within 250 pm depth only (for the used process condition) irrespective of energy partition 
due to EDM. 



Distance along x-axis (m) 
(a) 



(b) 

Fie 2J31 Effect of R«| on top surface temperature distribution in HSS workpiece 
rt „!««« T.9Ti« at v=1.0 mm Calonw X-X in Fie. 
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(b) 

Fig 2.32 Effect of R*. on temperature distribution along depth in HSS workpiece 
during EDDG process after 0.05 s at (a) x=1.0 mm and y=1.0 mm (Point A Fig. 2.21 b) 
(b) x=1.0 mm and y==0.5 mm (Point C Fig. 2.21b) 
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Time of machining (t) 

Effect of machining time on top surface temperature distribution is shown in Fig. 
2.J j(a) and 2.j3(b). A significant change in temperature is observed away from the spark 
radius. This is due to the fact that cut-off grinding temperature increases with the 
machining time. The peak temperature of about 1400 K is noted after machining time of 
0.04 s along X-X at y=1.0 mm (Fig. 2.21b). Temperature at spark center is increased by 
200 K after tour told increase in machining time. Fig. 2.34 shows the effect of machining 
time on depthwise temperature distribution. The effect is observed upto 500 ^im only. 



Distance along x-axis (m) 



(b) 


Fig 2J3 Effect of grinding time on top surface temperature dis^bution in HSS workpiece 
during EDDG process (a) along x-axis at y==1.0 mm (along X-X in Fig. 2.21b) (b) along y- 
axis at x“J.5 mm (along Y-Y in Fig 2.21b). 
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Table 2.1: Properties of workpiece used for Figures 2.5-111 (71) 


Workpiece 

AISI 1020 Steel 


19 mm X 6 mm X22 mm 

Wheel 

Resin bonded AI 2 O 3 

pw (kg/m^) 

7800 

C, (J/kg K) 

448.7 

k (W/mK) 

70.0 

To(K) 

298 


Table 12: Feed dependent values of and U, for AISI 1020 steel 
workpiece and resin bonded A1 1 0 j wheel during cut-off grinding [71] 


Vf (min/s) 

. _ Rvv. 

Uo (GPa) 

3.2 

0.31 

13.58 

4.0 

0.31 

13.58 

6.4 

0.29 

7.55 

9.35 

0.27 

6.2 

12.8 

0.23 

5.3 


Table 2.3: Properties of HSS workpiece used in present model [70] 


Size 

6 mm X 6 mm X 20 mm 

Wheel 

Metal bonded diamond wheel 

pw (kg/m^) 

8691 

p, (kg/m^) 

3520 

C, (J/kg K) 

419 

Cl (J/kg K) 

752 

Un (J/kgK) 

2735000.0 

Lv (J/kg K) 

15760000.0 

Tm (K) 

1965 

Tv (K) 

3100 

k(W/mK) 

40.2 

k,(W/mK) 

2000 

he (W/m^K) 

10000.0 (Kerosene) 

E (GN/m^) 

208 

Ok(/K) 

11.7 XIO^ 

aw(raV s) 

11 X 10 -^ 
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Table 2.4: Properties of workpiece used for Figures 2.15-2. 16 [70j 


Workpiece 

5 Cr Die Steel 

Pw (kg/m^) 

8593 

C,(J/kgK) 

425.0 

C, (J/kgK) 

752.0 

k (W/mK) 

48.5 

Lm(J/kg) 

2735000.0 

Lv(J/kg) 

15760000.0 

Tm(K) 

1965 

Tv(K) 

3100 


Table 2.5: Process parameters used for Figures 2.15-2.16 


Rvvs 

0.42 

1(A) 

6.5 

Ub (Volt) 

28 

R (pm) 

115 


Table 2.6: Process parameters used for Figures 2.17-2.34 


Vf (mm/s) 

10.29 

Rwg 

0.08 

Rws 

0.08 

Uo (GPa) 

10.55 

1(A) 

8 

Ub (Volt) 

40 

R (pm) 

125 

ton (ps) 

75 

toff(ps) 

150 
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2.5 Conclusions 

The present investigation shows that the temperature distribution in workpiece during 
EDDG can be simulated by means of FEM. This reduces the need for costly and time- 
consuming measurements to determine workpiece temperature. The proposed FEM 
approach takes care of the different length scales of the constituent process of EDDG 
(grinding and EDM) by using the superposition theorem. Thus, the 2-D temperature 
analysis of cut-off grinding and the axisymmetric analysis of EDM it is possible to 
predict the 3-D temperature distribution in the workpiece due to EDDG. Temperature 
studies based on 2-D temperature analysis of cut-off grinding and axisymmetric analysis 
of EDM compose favorably well with the experimental/theoretical results available in the 
literature. 

A sharp increase in temperature is observed around spark locations at the top surface 
of the workpiece. This shows that temperature-softening effects is limited to a thin layer 
around the spark flux. These are the well known facts in case of EDM but these facts 
were not verified in case ofEDDG quantitatively where sparking and grinding both are 
taking place simultaneously at the top surface of the workpiece. In EDDG the convection 
is taking place at the side surfaces due to presence of dielectric. I have found temperature 
and its softening effects quantitatively in case of EDDG where governing equations are 
different for EDM and grinding both with variable boundary conditions. The softening 
effect is significant upto 500 pm for the conditions taken. These results would be useful 
in the design ofEDDG operations. The model can predict temperature distribution in the 
workpiece including the effects of duty cycle, current, feed and time of machining if 
accurate values of energy partition due to both cut-off grinding and EDM, and for the 
conditions under consideration are known in advance. 



Chapter 3 

Thermal Stress Analysis of EDDG 


3.1 Introduction 

The concept of surface integrity has been widely accepted by researchers and 
industry experts in evaluating surface integrity of the machined surface. This term 
includes all aspects of the quality of surface, including surface finish, metallurgical 
damage and residual stresses [71], Machining operation produces a surface with a variety 
of surface alterations. The major types of alterations are; plastic deformation, micro-and 
macro cracking, untempered and overtempered martensite, phase transformation and 
residual stresses [77]. High temperature and/or high temperature gradients, and plastic 
deformation are the major cause of surface alterations. 

The thermal stresses in the workpiece are the best indicators of thermal influence 
on the surface integrity. Undesired residual stress distribution in a machined component 
results in lowering its dynamic strength, fatigue life, chemical resistance, magnetic 
properties and so forth. Need for theoretical study of thermal stresses arises because 
purely experimental investigation is too expensive and sometimes very difficult to 
perform. A theoretical method can take into account many principal effects 
simultaneously like mechanical, thermal, microstructural, and chemical. 

Heat generated within the machining zone causes a non-uniform temperature field 
in the workpiece leading to unknown material expansion of the workpiece [72]. Since 
bulk of the workpiece material under the top surface is still cool, expansion is constrained 
which produces compressive stresses in the sub-surface layer of the workpiece. The value 
of the thermal stress depends on temperature gradient and material properties. When the 
thermal stress reaches material yield stress value, permanent deformation occurs in the 
workpiece. As the workpiece cools, tensile stresses are created and remain in the 
workpiece. Therefore, a thorough examination of thermal stress development helps in 
understanding the relationship between machining conditions and thermal stress 


distribution. 
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Both the constituent processes of EDDG (cut-off grinding and EDM) generate 
heat, which can cause thermal stresses in the workpiece. The study of thermal stresses in 
the whole workpiece domain due to EDDG consists of thermal stresses due to both 
diamond cut-off grinding and EDM. 

In this chapter, an attempt is made to develop two independent FEM based 
models for the prediction of thermal stresses due to diamond cut-off grinding and EDM. 
Computation of thermal stresses is performed for various cut-off grinding and EDM 
parameters. 

3.2 Thermal Stresses in Cut-off Grinding 

Most of the thermal residual stress models available in literature [78-82] are 
limited to surface grinding for alumina wheels. Mishra and Prasad [78] applied finite 
element method to determine the residual stresses using the temperature field of Jaeger 
[35]. However, in their analysis the semi-infinite solid does not finally cool down to room 
temperature. Hence the use of the term residual stress is somewhat misleading. A 
comprehensive model for the quantitative prediction of thermal residual stresses in 
surface grinding is offered by Vansevenant [79]. An elastic-viscoplastic model developed 
by Li and Chen [80] has been used for predicting the residual stresses in surface grinding. 
But all these models assume that the line heat source or uniform rectangular heat source 
moves over the workpiece top surface. Mahdi and Zhang [81] have analyzed the thermal 
residual stresses induced by surface grinding with the aid of finite element method. 
Further, they have studied the coupled effect of thermo-plasticity and phase 
transformation on residual stresses [82], The model developed by Chen et al. [83] is 
capable of predicting the critical temperature as a function of grinding condition after 
which the tensile residual stresses are generated. 

In contrast to surface grinding, an analysis of thermal stresses during cut-off 
grinding is lacking. In this section, a numerical model is presented for the prediction of 
thermal stresses during diamond cut-off grinding. Finite-element codes have been 
developed to implement the thermal stress model for numerical computation. 
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3.2. 1 Mathematical Modeling 


The grain-workpiece interaction, which generates 
shown in Fig. 2.2. To make the analysis of thermal 


heat flux in cut-off grinding, is 
stress tractable a number of 


simplifying assumptions are required to be made as follows; 

Assumptions 


1 . A plane strain condition exists. 

2. The material of workpiece is homogeneous, isotropic and linearly elastic- 
perfectly-plastic. 

3. The yield point and thermal expansion coefficient are independent of 
temperature in contrast to real life situations, 

4. A moving heat flux of uniform rectangular profile is used to find the 
temperature distribution. 

5. No external forces act on the workpiece. Stresses occur only due to thermal 
expansion or contraction. 

6. Inertia and body force effects are neglected. 

7. The workpiece is assumed to be stress-fi-ee before grinding. 

Temperature model 

Temperature model discussed in section 2.2. (Fig. 2.3 and 2.4) is used to find the 
nodal temperatures in the workpiece domain. 

Thermal stress model 


The extreme temperature gradients that occur during cut-off grinding (in space 
and time) result in extreme nonuniformeties in the local thermal expansion of workpiece 
material (strain) which lead to high thermal stresses. The transient temperature 
distribution in the workpiece, obtained by solving the heat conduction equation (Eq. 2.2) 
along with the boundary and initial conditions (Eq. 2.3 and 2.4), is used as input for the 
calculation of thermal stresses. 

For thermal stresses due to cut-off grinding in a rectangular domain, the basic 
equations [76] and boundary conditions are as follows: 

Equilibrium equations 

The force equilibrium at a point leads to the following equilibrium equations 
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dx dz 


= 0 


d(J^ 

dx dz 


= 0 


(3.1) 


(3.2) 


Stress-strain-temperature relations 

The stress-strain-temperature relation for a linearly elastic material is 

{cr}=[i)]({.)-{™)) (3.3) 


where, the stress vector 


{cr}, the 


strain vector 


the vector (ot} and the elasticity 


matrix 



are given by: 



(3.4) 



(3.5) 



{ Eat AT 


(3.6) 



E 

(l + v)(l-2v) 


1- V V 
V 1 — V 
0 0 


0 

0 

l-2v 

2 


(3.7) 


Here, E is young’s modulus, v is Poisson’s ratio, at is coefficient of thermal expansion, 
and AT is temperature rise. 

Plastic zone 


To find the plastic zone, where the equivalent thermal stress is more than the material 
yield strength, the following expression for the deviatoric stress components are used. 
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^ cr. 
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^.X2 ~ ^XZ 


The equivalent (or effective) stress is given by 


CT, 


(s« + S Jy + )+ 2 s; 


The points in plastic zone are calculated by using the condition 
CTeq > CTy (plastic zone) 

Strain-displacement relations 

For small deformation, the strain-displacement relations are given by 
5u 5w 




( 5u 

dz 5x 


(3.8) 


(3.9) 


(3.10) 


(3.11) 


"‘ax’ 

Boundary conditions 

The domain of the problem with boundary conditions is shown in Fig. 3.1. Due to 
temperature effects, no rigid-body motion is allowed, but free expansion is possible. The 
bottom surface is not allowed to displace either in horizontal or in vertical direction. 


S2 


(tx = t,-0) 


1 SI (tx=tz = 0) 


(tx=tx = 0) 

S3 

S4 

►! 


(u * 0 and w * 0) 


Fig. 3.1 Boundary conditions applied to workpiece domain used for 
finding thermal stresses in cut-off grinding 
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Thus, 

u = 0 andw = 0atz = 0 (3.12) 

Traction on top and side surfaces is assumed to be zero i.e. 

t.= t. = 0 (3.13) 

The traction components tx and U are related to the components of the stress tensor by the 
following relation; 

t.=o-,jnj (3.14) 

where, nj are the components of the unit outward normal vector. 

3.2.2 Finite Element Formulation 


The finite element Galerkin formulation [73] is applied to convert the governing 
equations and boundary conditions to a set of algebraic equations. In this method, the 
displacements u and w are approximated by the some known functions. However, these 
approximations do not satisfy the equilibrium equations completely. The resulting error 
in them is made zero in an overall sense by using a set of fianctions Wx and Wj, called 
weighting functions. 

The resulting expression, called as the weighted residual integral, becomes 


J 



da. 


XZ 


dz 


r 






22 


8x dz 




dxdz=0 


(3,15) 


Here, the weight functions are chosen to have the value zero on the boundary on which u 
and w are specified. 

The use of divergence theorem leads to the following expression: 


\^Wf {a} dxdz = 0 p 

n 

Using the constitutive relation (equation 3. 3-3.7), this expression becomes 

n 

In the present analysis, 2-D eight-noded isoparametric element is used for 
descretization of the domain. Thus, displacement is approximated as 
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H=k]w 

where. 


(3.18) 


W ^ 

JVj 

0 

^3 

0 IV, 0 

^5 
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0 

0 fv, o' 

^ ^ L ^ 
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^2 

0 

0 N, 
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0 

0 

N, 0 N,_ 


is the shape function matrix and is the nodal displacement vector given by 

={wy Wj U2 W_2 Mj Wj U4 W4 Uj Wj Wg Uj Wy Ug Wg} 

Then the expression for strain becomes 


where, 

0 ^ 0 ^ 0 ^ 0 ^ 0 ^ 0 ^ 0 ^ 0 

kl- ? ^ 0 0 0 ft 0 0 f- 0 f- 0 ^ 

<w, av, dN, ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ 

IT dz Hx Oz dx dz dx dz dx dz dx Sz dxdzdx 


(3.20) 

(3.20) 

(3.21) 


For Galerkin method, we choose 


P;}=klH 

(3.22) 

where, {w' contains the nodal values of the weight Wx and Wz. Then 

{vir}=[r]^'*] 

where, are the nodal weight values. 

(3.23) 

Substitution of expressions (3.18-3.23) in equation (3.17) leads to 

the following finite 

element equations 


g ly‘J[sTK’]{sy=^'Z^‘JW’] 

(3.24) 

where, elemental coefficient matrix for stresses [STK‘ ]and elemental force vector for 

stresses {s 5 ^') are given by 

[S7X“]= 1 [af [^] [BjBxdz 

(3.25) 


k'}=L [bJ {m}dxcb 


(3.26) 
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Here, nem is the number of area elements and O' is the domain of typical area element. 
The integrals in equations (3.25-3.26) are evaluated by 3x3 Gauss quadrature integration 
technique. While calculating {m , the temperature at a gauss point is (Otained using the 
same interpolation functions as used for the displacements. 

The assembled finite element equations may be written as 

[GGRINDK]{GGRINDDU}= {GGRINDF] ( 3 . 27 ) 

where, [GGRINDK] is the global coefficient matrix, (GGRINDl j is the global nodal 
displacements vector and (GGRINDF) is the global right s.de vector for cut-off 
grinding. The final Eqs (3.27) are solved by Gauss-elimination method to find the nodal 
displacements after application of the displacement boundary conditions. 

Evaluation of thermal stresses 

As we have seen above, solution of the problem is obtained in the form of nodal 
displacements. From the nodal displacements, thermal stress at any point in an element is 
obtained by using Eq. (3.3). The equivalent stress is computed at each of the Gaussian 
integrating points and compared with the yield stress to identify whether yielding occurs 
or not (Eq. 3.10). This check is carried out for all the elements in the workpiece domain. 

3.2.3 Results and Discussion 

Present model has been applied to obtain the thermal stresses in HSS workpiece 

during cut-off gnndmg operation. The properties of HSS workpiece material are given in 
table 2.3. 

Fig 3.2 shows the variation of components of thermal stresses at top surface with 
grinding time due to uniformly distributed rectangular heat source. The nature of is 

compressive. The trend of can be explained by analyzing the temperature and 
temperature gradient of a point at the centre of top surface (x=0.003 m) of the workpiece 
for varying time. The value of is high for 0.1 s of grinding time because of high 
temperature (Fig. 3.3a) and also high temperature gradient with time fi-om initial time to 
0.1 s of grinding time (Fig. 3.3b). Further, the value of cr„ decreases significantly with 
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(a) (b) 


Fig. 3.3 (a) Top surface temperature variations with time (b) temperature variation of 
a point on top surface at x=0.003 m with time during cut-off grinding for U, =10.55 
GPa, Rwg =0.08 and Vf = 10.29 mm/s 

increase in grinding time from 0.1 s to 0.5 s because of decrease in temperature gradient 
within this time interval (Fig. 3.3b). After 0.5 s of grinding time, temperature gradient is 
not changing significantly and hence the value of is also not changing significantly 

(Fig 3.3b and Fig. 3.2a). The peak value of cr^ as 190 MPa after 0.01 s is decreased to 
120 MPa after 1.5 s. The values of and are very small as compared to cr^ (Fig. 
3.2b and 3.2c). 

Fig 3.4 shows the variation of thermal stresses along the depth with time. Here, 
the curves show that the non-uniformity in thermal stresses shifts from the top towards 
bottom with time of grinding. The trend of cr^ (Fig 3 .4a) can be explained by analyzing 

the temperature and temperature gradient with time of two points, one near the top 
surface (depth=0.0015 m) and other near the bottom surface (depth=0.019 m) along the 
depth at the centre of the workpiece (Fig. 3.5a and 3.5b). The temperature gradient with 
time near the bottom is more between 0.1 s to 0.5 s (Fig. 3.5a and 3.5b) hence high value 
of 3.4a) is obtained. Further, increase in grinding time does not increase the 

temperature gradient significantly and hence no further change in the value of cr^ is 
obtained. Near the top surface, temperature gradient is high till 0. 1 s of grinding time but 
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Fig. 3.4 Variation of (a) (b) cr„ (c) cr« with ttoe along depth 

at the centre of HSS workpiece during cut-off grinding for 
Vf-10.29 mm/s , R,,=0.08, U.= 10.55 GPa 
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Fig. 3.5 (a) Depthwise temperature variations with time (b) temperature 
variation of two points at different depths with time during cut-off grinding for 
x=0.003 m, Uo =10.55 GPa^ =0.08 and Vf = 10.29 mm/s 


its numerical value (approximately 250 MPa/s) is much less than the maximum value 
near the bottom surface (500 MPa/s, depth=0.019 m in Fig. 3.5b). The maximum values 
of components of thermal stresses developed either at the top surface or along the depth 
are less than the yield strength of workpiece material (3610 MPa). 

An important controllable parameter in cut-off grinding is feed given to the 
workpiece (Eq. 2.1) hence the effect of feed on variation of components of thermal 
stresses are analyzed as shown in fig 3. 6-3.7. 

Fig. 3.6 shows the variations of different components of thermal stresses with 
feed along top surface after grinding 5 ram depth of the workpiece. Here, it is observed 
that the nature of cr^ is compressive and varying with feed. Higher values of cr„ are 

obtained at higher feeds because of high temperature rise at top surface with higher feeds 
due to increase in qv,g (Eq. 2.1). The trends of cr^ and are also shown, but their 
quantitative values are very small (Fig. 3.6b and 3.6c) compared to . 

Fig. 3.7 shows the vanations of different components of thermal stresses with 
feed along the depth at the centre of the workpiece after grinding 5 ram depth. Here, it is 
observed that the thermal stress component cr^ varies with feed either near the top or 
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Fig, 3.6 Variation of (a) cr„ (b) and (c) era at top surface with Vr during 
cut-off grinding of HSS workpiece for R,g=0.08, 11,= 10.55 GPa, and 
machined depth 5 mm 
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Fig. 3.7 Variation of (a) cr„ (b) and (c) ct„ along depth at x=0,003 m with V, 
during cut-off grinding of HSS workpiece for R_-=0.fl8- IT = in «?«: np- 
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near the bottom surface (Fig 3.7a) Maximum increase of about 200 MPa in is noted 
near the bottom surface, by four-fold increase in feed (3.2 mm/s to 12.8 mm/s). Almost 
no change of cr^ and o ^ are noted with this change in feed (Fig 3.7b and 3.7c). 

3.3 Thermal Stresses in EDM 

High temperature generated due to high-density discharge thermal energy leads to 
not only thermal erosion but also formation of recast layer with micro-cracks on the 
machined surface. The presence of multi-layered heat affected zone and brittleness of the 
hardened layer have been reported [84-85] to reduce the fatigue strength of electro- 
discharge machined components. The computation of the temperature and thermal stress 
fields in the workpiece, therefore, is of considerable importance as far as surface integrity 
is concerned. 

Few available surface integrity models due to EDM are mainly based on 
metallurgical examinations of EDM’d surfaces using various experimental techniques 
such as optical metallography, scanning electron microscopy. X-ray diffraction [84, 86- 
88]. Lee et al. [87] have found the extent of surface damage with respect to more readily 
measured parameters such as pulse energy, and surface roughness. Kruth et al. [84] have 
studied the influence of workpiece material, electrode material and dielectric type on the 
composition and the metallographic phase of the white layer. However, there are very 
few theoretical approaches for determination of thermal stresses. 

Since the thermal action of EDM process affects the surface integrity of machined 
component, the analysis of temperature distribution is required for prediction of thermal 
stresses. Thermal stresses generated due to such a non-uniform temperature distribution 
can be analyzed for finding thermal residual stresses in the workpiece, which have direct 
relationship to the quality of the machined surfaces. 

In this section, a thermal stress model has been developed based on the finite 
element technique. Discharge current, pulse time, duty cycle, and spark radius are 
considered to be the main parameters affecting the thermal stresses. The principal aim o 
developing this model is to predict the nature of thermal stresses occurring during EDM. 
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3.3.1 Mathematical Modeling 

In EDM, two electrodes namely workpiece and tool, are dipped in liquid dielectric 
(kerosene) and physically separated by a gap, called inter-electrode gap. Due to random 
and complex nature of EDM process, the following assumptions are made to make the 
problem mathematically tractable. 

Assumptions 

1. The domain is considered as axisymmetric. 

2. The workpiece material is assumed to be homogeneous and isotropic. 

3. The material properties of the workpiece are temperature independent. 

4. The heat transfer to the workpiece is by conduction. 

5. Inertia and body force effects are negligible during stress development 

6. The workpiece material is elastic-perfectly plastic and yield stress in tension is 
same as that in compression. 

7. The workpiece is assumed as stress-free before electro-discharge machining 

8. Thermal stresses are evaluated only upto the time for which the transient 
temperature distribution is known. 

Temperature model 

Temperature model of EDM discussed in section 2.3 is used to find nodal temperature 
distribution within workpiece domain. 

Thermal stress model 

The extreme temperature gradients that occur during EDM (in space and time) 
result in extreme non-uniformeties in the local thermal expansion of workpiece material 
(strain) which lead to high thermal stresses. The transient temperature distribution in the 
workpiece, obtained by solving the heat conduction equation (Eq. 2.45) along with the 
boundary and initial conditions, is used as input for the calculation of thermal stresses. 

For the axisymmetric thermal stress analysis problem at hand, basic equations and 
boundary conditions are given below [76]. 

Equilibrium equations: 

The force equilibrium at a point leads to the following equilibrium equations 
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^ ^rr ^9B _ q 

dr dz r 


da^r 

dr dz 



(3.28) 


(3.29) 


where, Gnr, aee and Ozz are normal stresses and is shear stress. 

Strain-displacement relations: 


For small deformation, the strain-displacement relations in an axisymmetric domain are 
given by 

^ _du _ u _ dw _ du dw 

(3.30) 


dr 


-06 


dz ’ 


dz dr 


where Gn-, soe and Szz are normal strains and Srz is shear strain. Displacements along r 
and z-axis are u and w, respectively. 

Constitutive relations: 


The stress-strain-temperature can be written as: 

{o}=[D]{s)-{m) (3.31) 

The stress vector } , the strain vector } , the vector {m} and the elasticity matrix 
\p ] are given by: 

■" ^rr> ^09 > ^zz> ^rz 


v_^ 


] 


{w}: 


y| Eat AT 
l\ 


0 


l-2v 


[D]= 


(l + v)(l-2v) 


1-v V V 0 

V 1-v V 0 

V V 1-v 0 

0 0 0 


(3.32) 

(3.33) 

(3.34) 


(3.35) 
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Here, E is young’s modulus, v is Poisson's ratio, cu is coefficient of thermal expansion, 
and AT is temperature rise. 

Plastic zone: 

To find the plastic zone, where the equivalent thermal stress is more than the material 
yield strength, the following expression for the deviatoric stress components are used 

~ ^99 ~"z(p^rr ^ s! ^99 ) 

^ (3,36) 

S^=cr^--{cTrr+<r^+cx^) 

^rz ~^n 


The equivalent (or effective) stress is given by 


(3.37) 


where, a point in plastic zone is identified by the following inequality: 

cjeq > CTy (plastic zone) (3.38) 

Boundary conditions: 

The traction components tr and tz (mechanical stresses) are related to the components of 
the stress tensor by the following relation: 

U=(rijnj (3.39) 

where, nj are the components of the unit outward normal vector. 

Equations 3.28-3.29 are a set of elliptic equations in unknown displacements u 
and w, thus requiring boundary conditions along the boundary surface of the domain 
under consideration. The domain of the workpiece with boundary conditions is shown in 
Fig. 3.8. 

Displacements along and across surface AD are assumed zero. 

u = 0andw=0 on AD (3.40a) 

Since AB is an axis of symmetry, the normal component of displacement vector and 
shear component of traction vector are zero along the boundary. Thus, 
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u = 0 , = 0 on AB 


(3.40b) 



Fig. 3.8 Boundary conditions applied to workpiece domain used for finding thermal 
stresses in EDM 

Furthermore, since BC is a free surface, both the traction components are zero on this 
boundary. Finally, the boundary CD is considered sufficiently far away from the heat 
source so as to be stress free. Thus, the boundary conditions on BC and CD are: 
t = n t, = 0 on BC and CD (3-40c) 


3.3.2 Finite Element Formulation 

The finite element Galerkin formulation (73) is applied to convert the governing 
equations and boundary conditions to a set of algebraic equations. In this method, the 
displacements u and w are approximated by some known function. However, these 
approximations do not satisfy the equilibrium Eqs, (3.28-3,29). The resulting error m 
them is made zero in an overall sense by using a set of funaions W. and W, called as 

weighting functions. 
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Here, the weight functions are chosen to have the zero value on the boundary at which u 
and w are specified. 

The use of divergence theorem leads to the following expression: 

J {VWf {cr}2n: r dr dz = 0 (3.42) 

n 

Using the constitutive relation (3.3 1), we can write 


J ([r)]{s }“{'”}) 2 TT r dr dz = 0 


(3.43) 


n 


In the present analysis, 2-D eight noded isoparametric elements are used (Fig, 
2. 14). Thus, the displacement is to be approximated as 


where 

k]= 


N, 0 Nj 0 N, 0 N, 0 Nj 0 N, 0 N, 0 N, 0 

ON, 0 Nj O Nj 0 N, 0 N, O Nj 0 N, O N, 


is the shape function matrix and is the nodal displacement vector given by 

= {ui Wi U2 W2 Wj U4 W4 Wj Uq Uj W7 Ug W5) 

Then, the strain vector is given by 
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For Galerkin method, we choose 
fW. ' 


w. 


where, contains the nodal values of the weight and W^. Then 

{vfF}= [ b ] 


(3.44) 

(3.45) 


(3.46) 

(3.47) 

(3.48) 


(3.49) 


(3.50) 
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Substitution of expressions (3.44-3.50) in equation (3.43) leads to the following finite 
element equations 


f [w'f^TK-\[s\-=f_{w^y-} (35,) 

«-I 

where, elemental coefficient matrix for stresses ^TK ^]and elemental force vector for 
stresses W* are given by 

[sTr]=\{Bj[D\[B]rdrdz (3 5j, 

O' 

W’ }=jW {m}rdrdz (3 53) 

n 


Here, nem is the number of area elements and is the domain of typical area element. 
The integrals in equations (3.52) and (3.53) are evaluated by 3*3 Gauss quadrature 
integration technique. While calculating (m , the temperature at Gauss point is obtained 
using the same interpolation functions as used for the displacements. 

The assembled finite element equations may be written as 

[GEDMK]{GEDMU]= {GEDMF} (3.54) 

where, [GEDMK] is the global coefficient matrix, {GEDMU} is the global nodal 
displacements vector and {GEDMF} is the global right side vector. The final equations 
(3.54) are solved by Gauss-elimination method to find nodal displacements after 
application of boundary conditions. 

Evaluation of Thermal Stresses; 

As we have seen above, solution of the problem is obtained in the form of nodal 
displacements. From the nodal displacements, thermal stresses at any point in an element 
are obtained by using Eq. 3.31. The equivalent stress (Eq. 3.37) is computed at each of 
the Gaussian integrating points and compared with the yield stress to identify whether 
yielding occurs or not. This check is carried out for all the elements. 

3.3.3 Results and Discussion 

Present model has been applied to obtain thermal stresses in EDM operation. The 
properties of HSS workpiece material are given in table 2.3. The typical stress 
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distributions in EDM are presented first. Later, a parametric study of stress distributions 
is presented by varying the two parameters, i.e. duty factor and current. 

The thermal stress components are computed for each of the Gaussian integrating 
points within each element of the descretized domain (Fig. 2.14). All the Gauss points 
which are near the top surface of the considered domain and element size are found to be 
at 7 p.m below the top surface of the domain. Similarly, all the Gauss points which lie 
near the center line of the spark are found to be 7 pm away from the center line. Hence, 
typical planes, which are near to the spark, are chosen to show the results of stress 
distributions along the radius and along the depth for parametric study (Fig. 3.9 and Fig. 
3.12-3.15). 

Fig 3.9(a) shows the variation of thermal stresses along the radial distance at 7 pm 
below the top surface. It is noticed that a compressive stress region exists near the spark 
center (r = 0). The compressive nature of the stress components can be explained as 
follows. A part of the workpiece material under a spark expands more significantly due to 
higher temperature rise and causes compressive stresses because of the restraint from the 
material that expands less. Radial component (an-) is seen to be highest closely followed 
by tangential component {a-gg ). At the center, the values of each of these two components 
is 2750 MPa while that of Ozz is only 650 MPa. The shear component (CTrz ) ( 120 MPa) 
is much smaller compared to the other three components, and its sign changes a little 
away from the center of the spark. Almost no thermal stresses are observed beyond 400 
pm away from the center. 

Fig 3.9(b) shows the variation of thermal stresses along the depth. It is observed 
that the stress fields are tensile also in almost all cases. These curves also show that the 
undesired tensile thermal stresses are critical near the top surface. They are varying upto 
100 pm but gradually become zero beyond 400 pm. The more critical component is Ozz 
which varies from -675 MPa to +1300 MPa within 100 pm of depth. This shows that cr^ 
will cause very high thermal residual stresses near the top surface. 
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(a) 


(b) 


Fir. 3.9 Variation of stress components in HSS workpiece for Ub = 40 V, 1= 8 A, Rw 
=0.08, R=125 M-m, t„,=100 ^is, duty cycle=50%. Results after lOOps (a) along the radial 
distance at 7 pm below from the top surface, (b) along depth at r = 7 pm 


The zone where thermal stresses exceed the yield strength (Eq. 3.37) is 
determined after finding equivalent stress. The contours of equivalent thermal stress are 
shown in Fig. 3.10. The region starting fi-om the curve with Geq > 3610 MPa (near to 
3500 MPa curve in Fig. 3.10) moving towards the center of the spark would be under 
plastic condition where the stress will be greater than the yield strength of the workpiece 
(3610 MPa). But, the region near the spark where temperature is more than the melting 
temperature of the workpiece material (1935 K) (Fig. 3.11) is in molten state, and it is 
assumed that whole of the molten material gets ejected out. Therefore, no thermal 
stresses are developed in this region (Fig. 3.11) where workpiece temperature > 1935 K. 

Fig. 3.12 shows the effect of duty cycle on the thermal stress distribution along 
the radial distance. Details of on-time (U) and off-time (W) corresponding to each duty 
cycle is given in table 3.1. The curves clearly show that the radial and tangential stress 
components (cj-rr and age) increase with the increase in duty cycle value (Fig. 3.12a and 
3.12b). This is so because of the fact that heat supplied to the workpiece due to a spark is 
more with larger duty cycle and hence more temperature rise. Variation of radial stress 
component takes place upto 500 jmi (Fig. 3.12a) where as change in agg is only upto 175 
jam (Fig. 3.12b). Also, it is interestingly noted that shear stresses are changing their sign 
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Fig. 3.10 Equivalent isostress (N/m^) due to one spark within R/2 from the center 
(within one critical element) in HSS workpiece at Ub =* 40 V, 8 A, *0.08, R-125 
pm, ton="100 |is, duty cycle«50% 


Spark 



one 


Fig. 3.11 Isotherms (TK) due to oue spark within R/2 from the centre (within 

cntual element) in HSS workpiece at = 40 V, I- 8 A, R, -0.08, R-125 pm, 
t„,=100 ps, duty cycle=50%. Results after 100 ps. 
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within 200 |i.ni (Fig. 3. 12c). Almost no variation of is observed with duty cycle (Fig. 
3. 1 2d) Variation of stress components with duty cycle along the depth is shown in Fig. 
3.13. Here, it is noted that (Jrz and cTzz are not changing with increase in duty cycle, while 
other two are slightly changing. 



(a) W 



Fig. 3,12 Variation of stress components with duty cycle along depth in HSS workpiece at r 7 
(im from the center of the spark for Ub = 40 V , 1= 8 A, R* =0.08, R=125 pm, ton=100 ps (a) cs„ , 
(b) CTee , (c) CT„ , (d) Results after one on-time 

Table 3.1 Details of on-time, off-time and pulse time for different duty cycle 


Duty Cycle 

On-time 

i 

Off-time 

Pulse time 

((on-time/pulse time)* 100) 

(ps) 

(|4S) 

(H-j) 

200 

50 % 

100 

100 

60 % 1 

150 

100 

250 _ 

80 % 

400 

100 

500 







Stress (MPa) Stress (MPa) 
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Depth from top surface (m) 



Depth from top surface (m) 


(a) 


(b) 




Depth from top surface (m) 


Depth from top surface (m) 


(C) 


(d) 


Fig. 3.13 Variation of stress components with duty cycle along depth in HSS workpiece at r =7 
)im from the center of the spark for Ub = 40 V , 1= 8 A, R, =0.08, R=125 ^m, W=100 )is (a) ct,, 
(b) CTee (c) cr„ (d) <t^ Results after one spark on-time. 
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(a) (b) 



(c) 


(d) 


Fig. 3.14 Variation of components of tfaermal stresses with current along radius in HSS 
workpiece at 7 fjini below from the top surface for Uj, =40 V, I = 8A, R* =0.08, R=125 
^im, t«ff “100 ps, duty cycle = 50% (a) a „ , (b) Uee , (c) cr„ , (d) era. Results after 100 ps. 


The effect of current on the thermal stress distribution along the radial distance is 
shown in Fig. 3.14. When the current is high, the induced radial, tangential and transverse 
components of thermal stresses are also high and compressive in nature (Fig. 3.14a, 
3. 14b, and 3. 14d). In fact, the increase in temperature rise is directly proportional to the 





Stress(MPa) 


i08 


Thermal Stress Analysis of EDDG 




(a) 


(b) 




Depth from top surface (m) 


Depth from top surface (m) 


(C) 


(d) 


Fig. 3.15 Variation of components of of thermal stresses with current along depth in HSS workpiece 
at r = 7 pm for Ub =40 V, I = 8A, R, =0.08, R=125 pm, torr =100 ps, duty cycle = 50% (a) a,^ (b) cree 
(c) (j„ (d) CTn. Results after 100 ps. 

increase in current due to high heat flux value (Eq. 2.44) which leads to higher 
compressive thermal stresses. It is also observed fi-om Fig. 3.14c that the range of change 
of variation (with change in sign) of shear component of thermal stresses at the center of 
the spark is approximately in the same ratio as increase in the current. The variation of 
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these stress components with current along depth at r = 7 pm is shown in Fig. 3.15. Here, 
it is observed that the value ot any stress component increases in the same proportion as 
increase in current. The nature of variation of stress components is not changing with 
change in current value. 

3.4 Conclusions 

Two thermal stress modeis-one for diamond cut-off grinding and other for EDM, 
which are constituent processes of EDDG, have been developed to predict thermal stress 
fields in HSS workpiece by developing two finite element based codes. The results 
obtained serve to illuminate the development of damaging thermal stresses during 
diamond cut-off grinding and EDM. The thermal stresses developed during diamond cut- 
off grinding have been found to be below the elastic limit of the HSS workpiece material 
under the machining conditions used. Hence, no thermal residual stresses are expected 
during diamond cut-off grinding. It is also observed that due to a single spark substantial 
compressive and normal tensile stresses are developed in a thin layer. It is also found that 
the thermal stresses exceed the yield strength of the workpiece, which is mostly confined 
to an extremely thin zone near the spark. 



Chapter 4 

Thermal Analysis of Electro-Discharge Diamond 
Surface Grinding 


4.1 Introduction 

If EDDG is applied in surface grinding mode where whole workpiece and a part 
of the grinding wheel are completely dipped in dielectric, the process can be called as 
Electro-discharge diamond surface grinding (EDDSG) (Fig. 4.1). Similar to thermal 
analysis of EDDG as discussed in chapter 2 of this thesis, thermal analysis of EDDSG 
consists of thermal analysis of diamond surface grinding and EDM. Knowledge of 
thermal analysis of EDDSG with almost zero current is indispensable for analysis of 
EDDG in surface grinding mode. Metal bonded diamond wheels are normally used 
during EDDSG, which are structurally very different than the conventional wheels [71- 
72, 89-90], Hence, existing models pertaining to conventional wheels cannot be extended 
to diamond wheels. With this in view, a thermal model for diamond surface grinding is 
presented in this chapter. This model will be helpful in finding temperature distribution in 
the workpiece due to EDDSG. 

For grinding of steels and ceramics, embedded thermocouples and infrared 
detection technique have been extensively used [91] to measure workpiece temperatures. 
Beck [92] pointed out that the error in thermocouple reading due to the presence of 
thermocouple itself could be as large as fifty percent of temperature rise in certain cases 
when the thermocouple is placed normal to the heated surface as is often the case in 
grinding. The same difficulty persists in infi-ared detection technique as far as the 
problem of hole drilling in the workpiece for the purpose of placing the sensor is 
concerned. The temperatures measured by Kato and Fujii [93] at various depths from the 
ground surface were found to be of satisfactory accuracy and reproducibility. 

Most of the thermal models available in the literature are based on Jager’s moving 
heat source model [35-39]. Numerical methods have also been applied to find 
temperature distribution in the workpiece [49,51]. 
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Five parameters are important for determining temperature distribution; contact 
length (Ic), magmtude of heat flux going to the workpiece (qwg), distribution flmction ((|>g) 
for heat flux, energy partition (Rwg), and convective cooling (inside and outside the 
contact length) effect 


S2(hc 


Moving Heat Flux ( V w ) 



S4(T= To) 

Fig. 4.2 Thermal model used for EDDSG with negligible current 


Contact length: 

The effective contact length between the wheel and workpiece is one of the principal 
factors that determine the temperature distribution, deformation distribution as well as 
surface integrity of the workpiece. The magmtude of the contact length determines how 
the heat pulse is distributed over the workpiece surface. Many efforts [94-98] have been 
made to study contact length using experimental, analytical and numerical ^proaches. 
Gu and W^er [95] have performed microscopic study of contact length considering local 
deflection of grains, grain shape, chip formation, local heat transfer, surface roughness of 
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workpiece and so forth. Salje and Mohlen [96] have reported a relationship between 
geometric and kinematic contact length. Numerical methods have been utilized [97] for 
more accurate calculation of effective contact length. The predicted value of the real 
contact length vanes from one to ten times the geometrical contact length value [98], 
Snoeys et al. [32] have shown that the real contact length is also dependent on depth of 
cut as well as speed ratio of wheel and workpiece. 

The preceding studies show that at present there is no general formula that can 
predict the contact length applicable to all the grinding cases. Zhang et al. [94] have 
tabulated fourteen approximate formulae for contact length calculation. Most of them are 
based on experiments and their range of application is narrow. The relation between 
grinding conditions and the variation of contact length is complex. In the present model, 
geometrical contact length (Ic) is taken as equal to effective contact length for all 
computations, and it is given as follows 

l,^4Dd (4.1) 

The formula used for contact length by Eq. 4.1 ignores wheel-work deflection 
which will tend to increase wheel diameter but decrease depth of cut, sharpness of the 
wheel, coolant chemistry and its application methods, wheel and workpiece speeds as 
well as their properties. 

Heat flux profile and magnitude 

It is found that the spatial distribution (profile) of the heat source has a significant 
effect on temperature distribution in the workpiece. Usually, rectangular (uniform) heat 
flux distribution (Fig. 4.3) is assumed [35-38,51] for the thermal analysis of grinding but 
it is not supported by the temperature measurement. It is more appropriate to use a 
triangular heat flux distribution (Fig. 4.2) because the instantaneous material removal rate 
increases proportionally from the trailing to the leading edge of the grinding zone (Fig. 
4.1). Rowe et al. [99] compared the measured temperature distribution in the contact zone 
ivith theoretical one for three different profiles: rectangular heat flux distribution, 
triangular heat flux distribution, and square law heat flux distribution (Fig. 4.3). They 



Thermal Analysis o t Electro-Discharge Diamond Surface Grinding 


115 



Fig. 4.3 DifTercnt heat flux profiles used for thermal model of surface grinding (i) Rectangular 
heat flux, (ii) Right angle triangular heat flux, (iii) Square law heat flux distribution 


found that the best matching was with the square law heat flux distribution. Zhang and 
Mahdi [48] have shown that heat flux profile is neither rectangular nor parabolic nor 
square law. In fact, it can quite reasonably be approximated by a general triangular 
distribution with the apex located at a point between the inlet and the outlet of the contact 
bite (Fig. 4.2). Jen and Lavine [47] claimed that triangular heat flux results agree well 
with experimental results. In the present work, a general triangular heat flux profile (Fig. 
4.2) is used for diamond surface grinding. 

The magnitude of heat flux going to the workpiece during surface grinding is 
calculated by the following equation [72] 


f wg 




(^) 




(4.2) 


where, distribution function (j)g (Q for heat flux in surface grinding is given by: 
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1 

if ^<0 or I, 

if 0<^< + 

v2 J 

if ^ + <^<l, 

V 2 J 

Here, Rwg is energy partition due to grinding, Vw is table speed, is the distance of heat 
flux peak (A) from its centre (C) (Fig. 4.2), and Uo is specific grinding energy. 

Energy partitioning: 

Selection of energy partition (Rwg), which is defined as the ratio of energy carried 
by the workpiece to the total grinding energy, is of special importance in determining the 
thermal load on the workpiece. The value of this fraction is a controlling parameter in 
determination of the workpiece temperature rise and helps in understanding various 
thermal phenomena occurring within the grinding zone. 

Several researchers [99-111] have done extensive work to find the value of Rwg 
experimentally as well as theoretically during surface grinding. Ramanath and Shaw 

[100] studied the heat partition between an abrasive grain and the workpiece to obtain 

Rwg from the ratio of the geometric mean thermal properties (GMTP= -^kpc ) of the 

workpiece and the abrasive. The value of Rwg so obtained is 0.7 for aluminium oxide 
wheel and 0.51 for CBN wheel during pendulum and creep feed grinding of steel. Shaw 

[101] explained that the outstanding performance of CBN is not due to its high value of 
thermal conductivity alone but is really due to its relatively high value of GMTP, which 
in turn gives a low value of Rwg. 

Rowe et al. [102] determined energy partition by considering the grinding wheel 
(aluminium oxide) and workpiece (steel) as two sliders and assumed equal average 
temperature in the grinding contact zone. They have also shown that Rwg varies with table 
speed if table speed is more than 100 mm/s but the variation of Rwg with depth of cut is 
negligibly small for grinding of steel with CBN wheel. A substantially lower value of Rwg 
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was found for grinding with CBN and diamond wheels as compared to alumina and 
silicon carbide gnnding wheel with steel workpiece [103-105] because more grinding 

energy is carried away by the superabrasive grinding wheel due to its superior thermal 
properties. 

Kohh et al. [106] used the ‘trial and error’ method to obtain the value of R^g so 
that the values predicted by their model match with the experimental data on subsurface 
temperature during pendulum surface grinding of steel with aluminium oxide wheel. The 
value of Rwg so obtained has a very narrow variation (0.6-0, 7). Guo and Malkin [107] 
obtained the distribution of Rwg (earlier it was assumed constant) along the grinding zone 
under pendulum and creep-feed grinding by assuming the temperature on the workpiece 
surface equal to that on the composite surface at every point of the grinding zone. 

Energy partitions of 4 to 8.5 % for hardened bearing steel and 3 to 6 % for 
nodular cast iron are reported [108,109] under fine surface grinding with CBN wheel by 
using a temperature matching and inverse heat transfer method. These results indicate 
that vitrified bond CBN wheel may benefit from both high thermal conductivity of 
abrasive grain and cooling by grinding fluid due to porous structure of vitrified wheel. 
Cooling is improved because of the porous structure of the vitrified wheel enhances fluid 
flow. 

Kato and Fujii [110] estimated energy partition Rwg to the workpiece in surface 
grinding by measuring temperatures at various depths from the ground surface of the 
work-piece using PVD (physical vapor deposition) film method and top surface 
temperature by extrapolation. They reported that the value of the Rwg varies significantly 
from 0.3 to 0.8, depending on the combination of the workpiece material and the wheel 
material. For grinding of carbon steel, this was approximately 0.25 with the diamond 
wheel, 0.75 with the aluminum oxide wheel, and 0.45 with the silicon carbide and CBN 
wheels but independent of grinding parameters such as table speed (Vw) and depth of cut. 
For grinding with aluminum oxide wheels, the value of Rwg is about 0.5-0.6 for stainless 
steel and brass, and 0.7-0. 8 for cast iron and carbon steel. Ju et al. [Ill] have found that 
Rwg varies within contact length but this variation is negligibly small for fine grinding of 
steel with CBN wheel. Most of the suggested values are valid either for resin bonded or 
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vitrified bonded wheels. But information concerning the value of Rwg for grinding HSS 
with metal bonded diamond wheel is nonexistent in the literature. 

In the present work, the value of Rwg is chosen as 0.25 for computational purpose 
in diamond surface grinding which is close to the experimentally measured value in the 
reference [110]. 

Specific energy: 

A fundamental parameter derived from the power and grinding conditions is 
specific grinding energy, which is defined as the energy required to remove unit volume 
of material. It is obtained from the ratio of grinding power to material removal rate. The 
total specific grinding energy is considered to consist of chip formation, plowing, and 
sliding components. It is reported [71,72] that specific grinding energy is nearly identical 
for a diverse range of steels (low carbon steel to high speed tool steel) despite vast 
differences in workpiece hardness. 

Most of the temperature predictions available in literature are based on measured 
grinding forces from which specific grinding energy is determined for a particular 
grinding wheel, workpiece, and grinding conditions. It has been found that the specific 
grinding energy varies with depth of cut, table speed, work material, wheel material, and 
wheel conditions. It is very difficult to account variable specific energy purely from 
theory until unless nature of variation of Uo with each of these parameters is known in 
advance. It is also reported in literature [71,72,110] that the specific grinding energy lies 
between 60-85 GPa during grinding of HSS with CBN grinding wheel. 

In the present work, the value of Uo is chosen 55 GPa for fine diamond surface 
grinding which is close to the experimentally measured value in the reference [110] for 
surface grinding of steel with diamond wheel. 

Material thermal properties: 

In general, the material thermal properties such as specific heat capacity and 
thermal conductivity are temperature dependent. Most of the thermal models for 
grinding, available in literature, have assumed thermal properties independent of 
temperature. It is also reported that the values of thermal properties do not change 
sigmficantly with temperature and, hence, the average values can also be used. Using 
data of temperature dependent thermal conductivity and specific heat of HSS workpiece 
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[1 12]. the best-fit curves are obtained to find the nature of variation in polynomial 
equation form, which are used in the numerical calculations [Fig. 4.4]. These polynomial 
equations are given below; 

C = 286.97 + 0.578 T - 0.001 15 T" + 1 .5227' - 5.595 ;n0’'° J-* 

k = -39.296 + 0.5967 - 0.002287* + 4.6167' - 5.067** + 2.824 ZIO"'^ 7' - 6.244 XIO"'® 7^ 




Fig. 4.4 Variation of (a) thermal conductivity, (b) specific heat with temperature 
Dielectrics: 

The effect of surface cooling by grinding fluid was analyzed by Des Ruisseaux 
and Zerkele [38]. Cooling is not very effective for reducing local high temperature rise if 
extreme care is not taken for its application. It is critically important that the coolant is 
forced into the contact zone. Porous grinding wheels, such as electro-plated or vitrified 
bonded superabrasive wheels are very effective in reducing temperature in the contact 
zone. 

In the present work, dielectric acts as a coolant also. Generally, ho is a function of 
temperature, wheel speed and method of coolant application. No data is available in 
literature showing the nature of variation of ho with each of these parameters for grinding 
of HSS workpiece with diamond wheel. In this work, the value of convective heat 
transfer coefficient inside the contact zone is considered as only 60 % (0.6 he) of the 
outside value which is found after different computational experiments [Fig 4.8]. 
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4.3 Finite Element Formulation 


In the proposed model, the energy balance in the workpiece governs the heat 
transfer equations. The equations used for modeling and simulation include governing 
equations, boundary conditions, and initial conditions. 

Heat diffusion equation with temperature independent material properties is 
used to find temperature field due to grinding operation 


„ dT 5,.5T. d .,dT. 
pC, — = — (k — ) + — (k — ) 
dt dx dx dz dz 


(4.4) 


Here, T is temperature, p is density, k is thermal conductivity, and C, is specific heat of 


workpiece material. 

Energy transferred to the workpiece as heat input serves as thermal boundary 
condition on the top surface SI (Fig. 4.2). The bottom of the workpiece is assumed to be 
sufficiently remote from the top surface so as to remain at its initial temperature 
throughout the grinding operation. Heat loss to the coolant on the surfaces SI, S2, and S3 
is modeled using convective boundary condition. Thus, the boundary conditions are; 


q = -k'VT.n = -q^+h^T on Sl,S2,S3 


where. 


on 

on S2,S3 

T = Tf^ on S4 

The initial temperature can be taken as 
workpiece is completely dipped). Thus, 

T = Tf=T^ at t=0 

Galerkin Formulation 


'■ when t>0, (4.5) 

j 

temperature of the dielectric in which the 

(4.6) 


The Galerkin weighted residual method [73] is used to convert the governing 
differential equation (Eq. 4.4) into algebraic equations by using appropriate weighting 
functions. When the residue is made orthogonal to the weight function, the following 
integral equation is obtained. 
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Lkf 


at 


-.0 

.ox dx oz crjJ 


In the vector notation, this can be written as 
l^pC.^Wdxdz -\J^.{kVT)]Wctxdz =0 


(4 7) 


(4.8) 


The second term can be expanded using a vector identity. Thus, 

j^pC.^fVdxdz -Uv.{WkVT)]dxdz+j^7W.kVTdxdz =0 (4.9) 

The second term of the expression (4.9) is modified using the divergence theorem, which 
gives rise to certain boundary terms. Thus, the expression becomes: 

I, + l2 + l3=0 (4.10) 

where, 

I, =fpC, — Wdxdz (4.11) 

dt 


L= f VW.kVTdxdz =0 (4.12) 

JQ 

^3= ^ +lKJVTdr = 0 (4.13) 

Finite element approximation 

To obtain the finite element equation from expression (4.10), the domain is 
discretized into mm number of elements (Fig. 4.5). Further, over a typical element, the 
temperature is approximated using the unknown nodal values Ti and the known 
interpolation (or shape) functions N;* In the present work, 8-noded isoparametric 
quadrilateral elements are used. Thus, 

T{x,zj) = (4.14) 


where. 


{N-r={Nr 

N2 

N' .. 

n: / 

(4.15) 

and 

{r 

T' 

t; 

-T.‘] 

(4.16) 
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Fig. 4.5 The finite element mesh of the workpiece used in EDDSG 


In Galerkin’s method, the weighting functions are also expressed using the same shape 
functions. Therefore, 

nx.z) = {n-YH (4 17 ) 

where, the vector [W contains the nodal values of W which are known but arbitrary. 

Now, 


{vr}=[B-] {r-} 


(4.18) 




(4.19) 
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where. 
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(4.20) 


Further, 

?L 

dt 




where, i T * ^ contains the time derivative of the nodal value of T. 


(4.21) 


In the present work, the geometry is also approximated by the same shape 
functions. Thus, 






(4.22) 


where, the vectors (x'} and {z*]- contain the nodal values of x and z coordinates, 
respectively. 


To express the boundary terms in terms of the nodal quantities, approximation 
over the boundary is needed. Note that, as the domain is divided into nsm number of 
elements, the boundary automatically gets divided into nbe number of elements. Now, the 
approximation for W, T and qw over a t 5 q)ical boundary element can be written as. 





(4.23) 
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t = 'V3)<|r 


t: 


^2 ^3} 


^ ^ J 


(4 24) 


(4.25) 


Note that the vector {N*’} contains 1-D quadratic shape functions, which have a boundary 
approximation consistent with the 2-D approximation over the domain. 

Finite element equations 

Substitution of the approximation for temperature and the expression for weight 
functions (equations 4. 14-4.25) in integral (4. 10) leads to 

•=1 I J *«l 


where, the elemental capacitance matrix [C*] , the elemental conductivity matrix [A"] , 

the boundary conductivity matrix [aT*] , and boundary heat flux vector {/*) are given 
by 


[C']=J^, {n-Y dxdz 

W] = \^ k[B-Y [B-]dxdz 

{ivT^r 

{n‘’Y k‘kr 

4 


(4.27) 


Here, H is the domain of t)q3ical area element and Fb is the domain of typical boundary 
element (Fig. 4.5). The integrals in the expression (4.27) are computed using Gaussian 
quadrature with 3 points in each direction. The expression (4.26) when expressed in 
terms of global quantities takes the form 
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{Wl^lclfj} + {w}"[K]fr)= {Wf (F) (4.28> 

L J 

where, [C] is the global capacitance matrix, [K] is the global conduaivity matrix, PF] is 
the global heat flux vector, {T} is global temperature vector, and |r| is the time 
derivative of {T}. 

Since the nodal value of weight functions are arbitrary 

[C] jf| + W (4.29) 

Equation (4.29) represents a set of ordinary differential equations in the variable 
{T . These equations are solved after finite difference method (FDM as discussed in 
section (2.2.2). In diamond surface grinding, At is selected as time required to move 
forward exactly by one element length of the grinding domain in the direction of flux 
motion. 

4.4 Results and Discussion 

In order to assess accuracy and efficiency of the present finite element 
formulation, an example involving two different shapes of moving heat sources over the 
top surface of the workpiece are simulated and results are compared with the known 
theoretical results of Li and Chen [80] for conventional surface grinding. Table 4.1 gives 
the data related to thermal properties of grinding wheel and workpiece, and process 
parameters used in theoretical calculations of Li and Chen [80]. Due to limitation of 
computing resources, the domain size is chosen as one-fourth of the domain in reference 
[80]. Convergence study revealed that the mesh of 480 elements and 1577 nodes is 
adequate for this problem. Since no information is available about the actual shape of the 
heat flux, the results reported by the authors in this paper are obtained assuming two 
commonly used shapes: rectangular (uniform) and general triangular. The heat flux is 
distributed over 8 elements and moved forward exactly by one element after every step. 

A comparison of variation of surface temperature along the top surface with the 
results of reference [80] is shown in Fig. 4.6. The trend of the results matches well, but 
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there is discrepancy between the values. Tliis could be due to the difference in domain 
size and shape of the heat flux used in these two problems. 



Fig. 4.6 Variation of top surface temperature during surface grinding computed 
with heat source of different shapes and the results of Li and Chen {80J 


Table 4.1: Data related to thermal properties of grinding wheel and 
workpiece, and process parameters used by Li and Chen [SOJ 


Workpiece 

SAE1038 Steel 


(21 mm * 2.8 mm) 

Wheel 

AI 2 O 3 

pw (kg/m^) 

7880 

Cw(J/kgK) 

502 

k(W/inK) 

50,2 

(W/m-K) 

3000.0 (Water) 

V„ (nun/s) 

70.0 

V, (m/s) 

30.0 

D(nun) 

126 

d (mm) 

0.0625 

Uo (J/m^) 

68950 * 10* 

R*g 

0.70 

L (nun) 

2.8 

time (t) second 

0.2 
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Temperature distribution in the workpiece 

In this section, various results showing the effects of various process parameters 
on top surface temperature distribution in EDDSG with zero current are given. First, 
computational experiments are conducted to include the effect of temperature dependent 
properties of the workpiece on top surface temperature distribution in the workpiece but 
hardly any change in top surface temperature is observed (Fig.4.7). This shows that 
temperature generated for the set of conditions used is not sufficient to change the 
material properties significantly. Hence, for the prescribed machining conditions it is 
more reasonable to predict temperature based on average value of thermal properties. 
Process parameters used for the present calculation are given in Table 4.2. 


a 

a 

s. 

g. 

E 

O 

U 

m 

t: 

3 

(A 

Cl 

O 


550 

500 

450* 

400- 

350- 

300 ■ 


Based on average thermal properties 
-Based on temperature dependent thermal properties 




250 


"O.OOO 0.005 0.010 0.015 

Distance along x-axis (m) 


0.020 


Fig. 4.7 Effect of material thermal properties on top surface temperature at V* - 50 
mm/sec and d*=0.04 mm after 0.224 s 


Table 4.2: Process parameters used in the present model 


Workpiece 

HSS (21 mm * 2.8 mm) 

Grinding Wheel 

Metal bonded diamond wheel 

hc(W/m-K) 

10000.0 (Kerosene) 

Uo g/w?) 

55000 X 10^ 

Rwg 

0.25 
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Figure 4.8 shows temperature distribution at the workpiece top surface for 
different values of convective heat transfer coefficient inside the contact length. For 
simplicity, we first assume that inside convective heat transfer coefficient hi = f he, 
where y <I is a positive constant number (effective cooling factor). It is found that the 
critical value for HSS is 0.6. Above this value, there is no change in top surface 
temperature. It may differ from material to material. The effect is found to be 
insignificant because maximum temperature difference of 20-25 K is noted at the location 
where temperature peaks lie. 



Fig. 4.8 Effect of cooling on top surface temperature at V., = 50 mm/sec and 
d=0.04 mm after 0.224 s 

Figure 4.9a shows the top surface temperature distribution along the x-axis for 
different times. Here, At is selected as time required to move forward exactly by one 
element length of the grinding domain in the direction of flux motion. In Fig 4.9a the 
time interval for one element length movement is 0.007 s. After doing computational 
experiments it has been found that the peak temperature becomes constant after 0.056 s 
(after 8 element length movement of the heat source). It is also noted that the temperature 
gradient within the contact zone also becomes almost constant. The variation of 
temperature with time at a point (x=10.5 mm) on top surface of the workpiece is shown 
in Fig 4.9b. This shows the transient nature of variation of temperature over the 
workpiece top surface. 
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(a) (b) 

Fig. 4.9 Variation of temperature with time (a) at top surface (b) at a point 
on top surface ( x= 10.5 mm ) after different time for V* =“ 50 mm/sec and 
d=0.04 mm 


Figure 4,10 shows temperature distribution at the workpiece top surface for 
different depths of cut after grinding time t = 0.224 s when the table speed is constant. It 
is observed that the temperature at top surface is higher for larger depth of cut. From Eq. 
4.2 it can be seen that an increase in ‘d’ leads to higher heat flux at the contact zone 
which is the cause of increase in temperature at larger depth of cut. Here, it is observed 
that the temperature gradients which are the main cause of thermally induced residual 
stresses are increasing from 10.0 K/m to 85.0 K/m when d is increased from 0.01 mm to 
0.09 mm, respectively (between x=0.005 m to x=0.01 m). 



Fig. 4.10 Top surface temperature variation along x-axis for different depth 
of cut (d) at V w = 50 mm/sec after t = 0.224 s 
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Figure 4. 1 1 shows the effect of heat flux speed (table speed) for constant contact 
length. Here rise in surface temperature is noted when the length of heat source travel 
(Vw t) with different table speed is same (11.2 mm). The results show that the peak 
temperature of the top surface under the heat source increases with increase in table 
speed. This can be explained as follows. Equation (4.2) shows that the heat flux is 
proportional to the heat source speed. Hence by increasing table speed the heat flux will 
increase and hence the peak temperature will be more. The curves show that the 
temperature increases from 425 K to 600 K when the table speed is increased from 25 
mm/s to 100 mm/s. 



Fig* 4-11 Top surface temperature variation along x>axis for different 
table speed at d = 0.04 mm after table travel of 11.2 mm. 

Figure 4.12 shows the variation of top surface temperature for different diameter 
wheel (D) at a constant depth of cut. It is observed that the peak temperature is lower for 
higher diameter wheel. This trend can be explained as follows. With increase in grinding 
wheel diameter, the contact length increases which promotes more heat convection and 
decreases the amount of heat flux qwg (Eq. 4.2). 
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Fig. 4.12 Top surface temperature variation along x-axis for different diameter wheel at 
Vw = 50 mm/sec and d=0.04 mm after t * 0.224 s. 


4.5 Conclusions 

A FEM based numerical model has been developed to provide a theoretical 
determination of temperature distributions in HSS workpiece during EDDSG with zero 
current. This would help in predicting preliminary results, regarding temperature 
distribution in the workpiece, before going for actual costly and time-consuming 
experimental procedure for the measurement of workpiece temperature during actual 
EDDSG. It is found that the consideration of temperature dependent properties of the 
workpiece and non-uniform cooling consideration is not sigmficant. Top surfece 
temperature increases with increase in depth of cut and table speed. Bigger wheel is 
found to reduce the temperature rise at the workpiece top surface. 



Chapter 5 

Conclusions and Future Work 


5.1 Conclusions 

The present thesis work comprises finite element analysis of EDDG process for 
finding temperature distribution in the workpiece only. Temperature distribution in the 
workpiece due to EDDG has been obtained by thermal analysis of diamond cut-off 
grinding and EDM. Thermal stresses have also been obtained for each constituent process 
of EDDG such as diamond cut-off grinding and EDM. EDDSG has been analyzed for 
temperature distribution in the workpiece when current is switched oflf. Numerical 
experiments have been conducted to study the influence of process parameters on the 
process performance. On the basis of these computational investigations, following 
conclusions can be drawn: 

❖ Thermal Analysis of EDDG 

■ 3-D temperature distribution in the workpiece during EDDG can be obtained by 
superposing results of the 2-D thermal analyses of diamond cut-off grinding and 
axisymmetric thermal analysis of EDM. 

■ Temperature fields of 2-D thermal analysis of cut-off grinding and axisymmetric 
thermal analysis of EDM compare favorably well with the experimental/theoretical 
results available in the literature. 

■ The thermal model for diamond cut-off grinding predicts the temperature 
distribution in the workpiece fairly accurately including the effects of feed, 
convective heat transfer coefficients of dielectrics and thermal properties of 
diamond grinding wheel if accurate values of specific energy and energy partition 
for the process conditions are available. 

■ The thermal analysis of EDM predicts the temperature distribution in the workpiece 
fairly accurately including the effects of on-/oflf-time, phase change, duty cycle if 
accurate values of the energy partition and spark radius for the process conditions 


are known. 
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■ A sharp increase in temperature is obser.-ed around spark locations after 
superposition of grinding and EDM temperature distributions at the top surface of 
the workpiece. 

“ Temperature gradients due to EDDG are found to be very high in a. thin surface 
layer adjacent to spark locations. This shows that temperature-softening effects are 
limited to a thin layer around the spark flux. 

■ Effects of feed and current on temperature distribution during EDDG are 
significant but no change is observed by varying duty cycle. 

♦> Thermal Stress Analysis of EDDG 

■ Two thermal stress models developed, one for cut-off grinding and other for EDM, 
are capable to determine the thermal stresses due to each constituent process of 
EDDG. 

■ The thermal stresses developed during cut-off grinding have been found to be 
below elastic limit of the HSS workpiece material for the conditions used. Hence, 
no thermal residual stresses are found during diamond cut-off grinding under the 
machining conditions used during the reported work. 

■ The thermal effect of spark is limited to a localized area around the spark at the top 
surface of the workpiece. Hence, thermal stresses have been found in the plastic 
range within a small thickness around the crater created by the spark. 

Thermal Analysis of Electro-Discharge Diamond Surface Grinding 

■ A FEM based numerical model developed for simulation of temperature 
distribution in the workpiece during EDDSG with negligible current (case of 
diamond surface grinding in which workpiece is dipped in dielectric) would reduce 
the need for costly and time-consuming experimental procedure for the 
measurement of workpiece temperature during grinding. 

■ The results presented would be useful for prediction of temperature distribution in 
the workpiece in combined process of diamond surface grinding and EDM called 
as EDDSG. 

■ The effect of cooling on temperature rise in the workpiece is found to be less 
significant during diamond surface grinding. 



Conclusions and Future Work 


135 


5.2 Scope for Future Work 

In the present study, an attempt has been made to determine 3-D temperature 
distribution in the workpiece during a new hybrid machining process called Electro 
discharge diamond grinding (EDDG). In this thesis, thermal aspects relating to temperature 
distribution and thermal stresses have been studied for deep understanding about EDDG 
process. Some of the issues, which can be examined in future, are summarized as follows: 

1. A comprehensive 3-D thermal model is needed to find temperature distribution in 
the workpiece during EDDG in cut-off grinding mode considering random 
generation of sparks at the top surface. 

2. A fast response temperature-measuring instrument and a suitable dynamometer 
need to be employed to measure temperature and cutting forces respectively. The 
temperature dependent specific energy can be estimated using measured cutting 
forces and then actual strength of heat flux supplied to the workpiece can be 
determined. Again, energy partition to the workpiece during EDDG can be 
obtained by temperature-matching methodology using measured temperatures at 
various depths. 

3. An integrated theoretical approach could be tried for thermal analysis of EDDG 
process where the effect of combined phenomenon of abrasion and EDM on 
temperature distribution in the wheel and workpiece can be treated together. First, 
the energy partitions into wheel, workpiece, and dielectric are to be found and then 
using these data, temperature distribution in the wheel and workpiece can be 
determined. 

4. The feasibility of EDDSG could be explored experimentally to generate flat 
workpiece surfaces. 

5. In addition to temperatures, a rigorous theoretical analysis, as well as experimental 
measurement of residual stresses produced in the workpiece are completely 
lacking. Hence, a surface integrity model needs to be developed which can predict 
transient temperature field, residual stresses and thermal damage in the workpiece. 

6. Because of the complexity of EDDG process, the optimization of process 
parameters as well as optimal control of the process will be of considerable interest. 
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Appendix A 

Derivation for heat flux due to Gaussian distributed 
spark heat source 





Considering the effect of spark upto 6-a limit i.e. 
dcr = 2R or CT = R/3 

The heat flux qws (r) at a radius (r) can be calculated if maximum intensity qo at the axis 
of a spark and its radius (R) is known. 

^v,,('') = ^oexp(-4.5 AI.l 

At ^=0 , ri=qo and at ^=R , Ti=qo e"^’’ 

Now from equation Al. 1 

In— = -4.5^ or =- — ln— A1.2 

If it is assumed that total power of each pulse is to be used only by one spark then 
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Total volume of the of the spark is given by: 

V=tw Ub 1= VI + V 2 A1.3 


Where 

Volume after rotating OOfC D about 0 02 fVj ) = z e~^'^ = 0. OllrdV AI .4 


Volume after rotating O^ O 2 C about OiOlfV^ ) = tc 



= ^J . 

Jq,e 


f n2 

Ir,” 


= 0.20867tR^q^ 


dq 


A1.5 


Putting Al. 4 and A1.5 in A1.3 

fwUbI = O.OIIkR^ q^^ 0.2086nR^q^ 

OR 

^ _4.45R^,U, I 

TT 

Putting the value of qo from equation A 1.7 into equation Al.I 

TV — V 1 

Where Ub is breakdown (discharge) voltage and I is current 


AI.6 


A1.7 


A1.8 
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